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Abstract - The spallation residues produced in the bombardment of 5e Fe at 1.5, 1.0, 0.75, 0.5 and 0.3 A GeV 
on a liquid-hydrogen target have been measured using the reverse kinematics technique and the Fragment 
Separator at GSI (Darmstadt). This technique has permitted the full identification in charge and mass of 
all isotopes produced with cross-sections larger than 10~ 2 mb down to Z = 8. Their individual production 
cross-sections and recoil velocities at the five energies are presented. Production cross-sections are compared 
to previously existing data and to empirical parametric formulas, often used in cosmic-ray astrophysics. 
The experimental data are also extensively compared to different combinations of intra-nuclear cascade and 
de-excitation models. It is shown that the yields of the lightest isotopes cannot be accounted for by standard 
evaporation models. The GEMINI model, which includes an asymmetric fission decay mode, gives an overall 
good agreement with the data. These experimental data can be directly used for the estimation of composition 
modifications and damages in materials containing iron in spallation sources. They are also useful for improving 
high precision cosmic-ray measurements. 



I. INTRODUCTION 



The spallation cross-sections of nuclides such as Fe have been historically studied to understand the propa- 
gation of cosmic-ray ions in the Galaxy, and to determine the composition of the Galactic Cosmic Ray (GCR) 
source [l|- @. Galactic cosmic rays constitute a superthermal gas that is partially confined in the Galaxy by 
interstellar magnetic fields with some leakage into the intergalactic medium. While propagating in the Galaxy, 
cosmic rays pass through the interstellar medium and some primary cosmic ray nuclei spallate into secondary 
cosmic ray nuclei. As measured by instruments in the solar system, the composition includes both primary 
cosmic rays whose abundance is depleted by spallation, and secondary cosmic rays produced by spallation. As 
a result of spallation during propagation, certain elements in the GCRs are far more abundant (often by orders 
of magnitude) than in solar system material. Examples of these "secondary elements" include Li, Be, B which 
are mainly spallation products of C and O, and Sc, Ti, V and Cr which are mainly spallation products of 
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Fe. Conversely, those elements where the abundance of heavier elements is much smaller, and hence have very 
small secondary contributions are "primary elements." Prominent examples include C and and Fe. Provided 
the spallation cross-sections are known, the abundance of secondary elements relative to primary elements are 
a measure of the amount of material cosmic rays traverse in the Galaxy. This in turn constrains astrophys- 
ical models of cosmic rays in the Galaxy. It is possible to correct abundance measurements for propagation 
back to the "source," that is, to determine the composition of the material that became the cosmic rays. The 
secondary-to-primary ratios combined with the cross-sections determine the amount of material traversed dur- 
ing propagation in the Galaxy; the amount of material traversed, again with the cross-sections, is then used to 
correct the measured abundances to the source abundances. Thus, uncertainties in the cross-sections are more 
significant than any details of the astrophysical models. (The exception to this generalization are the unstable 
secondaries.) In recent years, new high resolution elemental and isotopic measurements have become available 
(i.e. the ACE [Io| and Ulysses [Tlj space missions), including measurements in the iron region. The main source 
of uncertainties in determining both cosmic-ray secondary production and source composition using these data 
are uncertainties in the nuclear cross-sections. The interstellar medium is composed ~ 90% by number of H 
atoms and ions. Most high resolution measurements are of cosmic rays with energies per nucleon in the inter- 
stellar medium of ~ 0.5 to ~ 1.5 GeV. The cross sections reported here are thus directly applicable to improved 
interpretation of high-precision cosmic-ray measurements. 

Spallation reactions have also gained a renewed interest with the recent projects of spallation neutron sources 
and accelerator-driven sub-critical reactors systems considered for the transmutation of nuclear waste (Acceler- 
ator Driven Systems (ADS)). In these systems, a high-intensity proton beam of energy around 1 GeV is guided 
on a spallation target made of a high-mass material. In ADS, neutrons produced in the spallation target are 
used to maintain the reactivity in the sub-critical reactor where nuclear waste can be transmuted. The proton 
beam under vacuum in the accelerator has generally to cross a window before entering the spallation target. 
As it is continuously submitted to the proton beam irradiation, it is one of the most sensitive parts in ADS or 
spallation-neutron-source design. Among the problems created by the proton irradiation are the changes in the 
chemical composition of the window material and embrittlement created by gas production and atomic displace- 
ments (DPA) in the crystal lattice. A large range of materials have been studied for this window and, in most 
of the projects, martensitic steels composed at 90 % of Iron (with also substantial quantities of Chromium and 
MoUbdenium) have been retained due to their resistance to thermal constraints and radiation effects. Therefore, 
it is important to have a good knowledge of the production cross-sections of spallation residues in Iron and of 
their recoil velocity. 

In recent years, an important effort has been undertaken, mainly under the framework of the HINDAS 
European proje ct fl2l| , to collect a comprehensive set of high-quality spallation data regarding the production 
of neutrons [lj, Il4|. light charged particles [l5| and residual nuclei. The general goal is to better understand 
the reaction mechanisms in order to improve the models that are then implemented into high-energy transport 
codes. These codes, validated on experimental data, can afterwards be used to reliably predict all quantities 
needed for the design of ADS or spallation sources as neutron production, activation or damages. 

As concerns residue production, up to now the emphasis was put on spallation reactions on heavy nuclei. 
Isotopic cross-sections of residues produced in the reactions 197 Au + P_H, G3 at 800 A MeV , 208 Pb + p at 1 
A GeV and 500 A MeV [H d M, EI, 238 U + p at 1 A GeV [H, \2j$T 238 U + d at 1 A GeV 0, Hg have 
already been measured using the reverse-kinematics method at GSI (Darmstadt). In this paper we present 
new experimental results concerning the isotopic production cross-sections and recoil velocities of spallation 
residues in the reaction 56 Fe + p for five energies of the iron beam (0.3, 0.5, 0.75, 1.0 and 1.5 A GeV). This 
measurement is the first consistent set of data on isotopically identified residues on a large energy domain and 
for a light nucleus of practical interest. The comparison of the obtained data with various models, some of them 
being quite successful for heavy systems, allows testing their predicting capabilities for light nuclei and their 
dependence on beam energy. 



II. EXPERIMENTAL METHOD 



A. Experimental set-up 

In October 2000, an experiment was performed using the reverse kinematics at GSI in Darmstadt, Germany. 
A primary beam of 5e Fe was delivered by the heavy-ion synchrotron SIS at energies of 0.3, 0.5, 0.75, 1.0 and 
1.5 A GeV and directed onto a liquid- hydrogen target designed and built in the Laboratoire National Saturne 
(Saclay, France) 

The liquid- hydrogen thickness was &7.3mg/cm 2 contained by titanium windows of 20pm each. Two additional 
Ti foils were used to isolate the vacuum around the target from the vacuum of the beam pipe for security reasons 
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so that a total of 36mg/cm 2 of Ti contributes to the empty-target counting. Measurements were repeated with 
an identical empty target in order to subtract the production on the titanium container from the measured 
yields of residual nuclei. The contribution of these walls to the counting rates was below 10 % for the main part 
of the residues and below 20 % for the lightest ones. 

The time structure of the primary beam was a pulse of 6 s every 12 s, and the intensity was limited to 10 7 
part/spill. This beam intensity was measured using a secondary-electron emission monitor (SEETRAM) [27j 
calibrated at the beginning and at the end of each set of measurements at a given beam energy. This was done 
at low counting rates with a plastic scintillator as absolute reference. 

Residual nuclei produced in the reaction with the target were focused in the beam direction and analyzed with 
the FRS (Fragment Separator) [28| operated as an achromatic magnetic spectrometer. Fig. [T] is a schematic 
diagram of the experimental setup showing the four large dipole magnets and the essential detector equipment. 



MWPC 



Scintillator 



SEETRAM 



Dipole magnets 



Liquid H2 
target 

SO 




ToF (37 m),- 




» V 



Scintillator 



MUSIC Focal P lane 
ionisation S4 
chamber 



Dispersive mid-plane 
S2 



FIG. 1: Schematic layout of the FRS fragment spectrometer. Fragments are analyzed by the four large dipole magnets. 
Scintillators at S2 and S4 measure the time of flight over the second half of the spectrometer as well as the horizontal 
positions in the dispersive focal planes at S2 and at S4. The MUSIC detector (ionization chamber) gives information 
about the energy loss of the fragment. Multi-Wire Proportional Chambers (MWPC) are used for beam tuning and 
removed for production measurements. 



Due to their relativistic energies, the fragments produced in this experiment are fully stripped. The horizontal 
positions of these ions and a time of flight (ToF) were measured with two plastic scintillators, one located in 
the intermediate dispersive plane S2 and the other one installed at the final achromatic focal plane S4. The 
signal from the scintillator at S4 was used as the trigger for the acquisition of all detectors. The nuclear charge 
Z was determined using a multiple-sampling ionization chamber (MUSIC). The energy loss in the gas produces 

a signal proportional to 4s-, allowing the determination of Z with a resolution of AZ = 0.3 (FWHM) charge 
units. 

The knowledge of the horizontal positions of the ions determines precisely the radii pi and p2 of their 
trajectories in the two magnetic sections of the spectrometer. An absolute calibration is obtained with the iron 
beam detected in specific measurements at low intensity. Together with the magnetic field strengths in the 
dipoles measured with Hall-cffcct probes, the magnetic rigidities Bpl and Bp2 can be determined for each ion. 
Therefore, a total identification of the nature of the ions could be performed from the relation : 



A _ eBp 

Z m u c/3j 

where m u was the atomic mass unit and [3 7 were deduced from the experimental time of flight. Note that 
in this formula we have replaced the mass of the {A, Z) ion by A.m u which means neglecting binding energies 
compare to nucleon masses. 

The FRagment Separator has a momentum acceptance of ±1.5%. Therefore, about 18-20 different settings 
of the FRS were needed to cover the complete velocity distribution of all the ions. Figure [2] shows the complete 
fragment coverage in the Z vs. A/Z plane for 1 GeV per nucleon 56 Fe on the hydrogen target. The plot was 
made by adding histograms from individual settings, each one normalized to the dose of the primary beam. 
Fragments are well resolved and easily identifiable in this plot down to lithium. However, for the lightest 
elements the transmission of the spectrometer is very low, necessiting a dedicated method of analysis. This has 
been done only at 1 GeV per nucleon and reported in a separate paper J2f|. Therefore, we show in this paper 
results of the production cross section and recoil velocity only down to Z=8-10, depending on the beam energy 
considered. 
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FIG. 2: Complete isotope coverage in Z vs. A/Q (actually identical to A/Z) for 1 A GeV 56 Fe on the liquid-hydrogen 
target. The plot is built from data of overlapping settings, normalized to the primary beam intensity and superimposed. 

B. Data analysis 

The fragments are first identified in Z using the ionization chamber, taking into account the position and 
velocity dependence of the energy-loss signal. The velocity distribution of the fragments is obtained with high 
precision using the time-of-flight and magnetic-rigidity measurements. The experimental time-of-flight between 
the intermediate and the final focal plane is precise enough for an unambiguous identification of the fragment 
mass. After identification of the isotope, a more accurate value of the longitudinal velocity can be deduced from 
the magnetic rigidity in the first part of the spectrometer using relation 1 . 

Assuming that the reaction takes place at the center of the target, the fragment velocity is corrected for the 
energy loss in the target and transformed into the reference frame of the projectile at rest. A measurement of 
the recoil velocity of the fragments is thus obtained in that frame. To obtain the production cross-section of a 
given isotope, it is necessary to reconstruct the full velocity distribution by adding the partial ones measured 
in different settings, with the proper normalization. An example of the velocity distribution for 38 K is shown 
in Fig. [31 For this isotope, five different settings of the FRS were needed in order to reconstruct the complete 
velocity distribution. 



Due to potential damages in the detectors, isotopes having a magnetic rigidity too close to the beam one could 
not be measured. This is why the detection of 54 Mn was not possible in this experiment. For the same reason, 
some settings of other isotopes could not be obtained, leading to truncated measured velocity distributions. In 
that case a fit by a Gaussian function excluding the truncated zones was used to reconstruct the full distribution 
and then determine the total cross-section, the mean value of the velocity and its variance. In the case of a 
truncated zone in the velocity distribution too large to have a converging fit, the parameters of the Gaussian 
were constrained using the neighboring isotopes. The reconstruction procedure leads to an uncertainty on both 
the velocity determination and the isotope production cross-section. These uncertainties have been estimated 
by taking into account the fluctuation of the reaction point in the target and by doing reasonable variations of 
the fitted parameters for several groups of isotopes. 



5 



0.35 



0.3 - 




FIG. 3: Longitudinal velocity distribution of K detected as a residual nucleus at 1 A GeV and expressed in the rest 
frame of the iron beam. Five different settings of the FRS were needed to reconstruct the complete distribution. The 
yield (here in arbitrary units) detected in each setting is normalised by the number of incident iron nuclei and corrected 
for the acquisition dead time. 



C. Corrections and uncertainties 



The isotopic production cross-section of each spallation residue <j{Z, A) was obtained from the difference 
between the yield measured with the hydrogen target (Yh{Z, A)) and the yield measured with the empty target 
(Y e (Z, A)), each of them corrected for their dead time (correction factor / T ^(/ Te )) and normalized to the number 
of incident iron nuclei Np e h{Nfc e)- 



{Z,A) ■ f T H Y e (Z, A) ■ f Te \ / e • f trans ' f sec , s 

a(Z ' A) = {—^Th ) ^ (2) 

The cross section is finally obtained after a division by the number of hydrogen nuclei per surface unit Njj 
and with additional corrections due to the detection efficiency (/ e ), the transmission of the FRS (ftrans) and 
the secondary reactions (/ sec ) estimated for hydrogen events. It was determined that, even at the lowest energy, 
a correction for possible change of charge state is not necessary. 

Losses of events due to the dead time of the experiment, mainly due to the acquisition capability, are estimated 
for each run from the ratio between the free triggers measured on a scaler of high counting-rate capability and 
the number of recorded events (or accepted triggers). During the experiment, the counting-rate conditions were 
kept so that this correction never exceeded 30%, and was most frequently smaller for detection at magnetic 
rigidities substantially different from the beam rigidity. 

An estimation of the global detection efficiency f e including the detailed analysis of all needed information 
can be obtained from the difference between the number of accepted triggers and the final number of events 
that have been analyzed. An event can be analyzed if all the elements required have been registered without 
any problem: position at the two focal planes, time of flight and energy loss in the MUSIC detector. The trigger 
signal obtained by a narrow coincidence on high signals produced by highly ionizing particles is here supposed 
to identify a true heavy ion with a probability of nearly 100%. In almost all settings this efficiency was in the 
range 96-99%. 



() 



Corrections due to secondary reactions in the target and in the layers of matter on the trajectory of the frag- 
ments (mainly the plastic scintillator of 3 mm thickness at S2) were calculated following the method described 
in [3(| as previously used in other similar experiments [l^ . If a second reaction occurs in the target, the initially 
produced ion becomes lighter, so that cross sections of light ions are artificially increased (and the one for the 
corresponding heavy ion decreased). If a reaction occurs in the plastic at S2, the spallation ion will most often 
be out of the narrow magnetic rigidity acceptance in the second part of FRS and so will be lost at S4. Total 
nuclear interaction cross sections for the different fragments were estimated using the parametric formula of 
Kox et al [3l|. The maximum value (8%) of this correction factor is obtained for the secondary reactions in 
the target leading to the lightest evaporation residues. It decrease to zero for heavy residues. The correction 
due to the lost in the scintillator if a reaction occures is of the order of 3.5% and was taken into account (as a 
function of the nature of the ion and of it's mean energy). The attenuation of the beam flux inside the finite 
target thickness was also taken into account in this correction and is equal to -2% for a reaction cross section 
of 700 mb. 



The transmission correction is the most important factor concerning losses in the detection. Due to its 
geometrical characteristics and the ion optics, the FRS has only an angular acceptance of 15 mrad around the 
beam axis, and a large number of the fragments analyzed in this experiment have an angular distribution at the 
entrance of the FRS larger than this acceptance. An evaluation of the fraction of the residual yield not detected 
in the experiment had to be made from the measured velocity distribution of the fragment as it is described 
in [32| . Considering that, in the projectile reference frame, the emission of the fragments can be described as 
a 3-D Gaussian distribution around a mean longitudinal recoil, the width of the angular distribution in the 
laboratory frame can be obtained from the longitudinal velocity distribution measured in the experiment: 



« ^4 (3) 

where < uii > is the mean value and ct(i>||) the width of this distribution for evaporation residues of a given 
mass. 

The transmission through the FRS can be parameterized as : 



where a e ff(x2, x4) is the effective angular acceptance of the FRS as a function of the ion positions x2 and 
x4 respectively at the intermediate S2 and the final focal planes S4. This angle was calculated with the code 
described in reference [32| using 15 mrad as the maximum angular acceptance when the ion optics is the most 
appropriate. 



The transmission factor varies from 1 (no correction) to 0.4 for the lightest fragments that have a much 
larger angular distribution (see Fig. [4]) for the three highest energies. Various reasonable assumptions on the 
calculation of a e ff(x2,x4) lead to uncertainty estimations on T of 1% to 15% for the lightest evaporation 
residues. However, the analysis has revealed that at 500 and 300 MeV/A, the magnetic optics settings used 
during the experiment was not optimal and that the maximum acceptance of the FRS was reduced to 9.15 
mrad. This value has been taken into account in the transmission factor leading to much larger corrections for 
these two energies as it can be seen in Fig. 01 

For the absolute normalization, the precision on the target thickness has been studied in previous experiments 
[33| and is estimated to be 2.5%. The absolute numbers of incident ions Np e h and Np e e for runs with the 
hydrogen target and the empty target respectively are obtained from the SEETRAM calibration with an absolute 
error estimated to be 2.8%. 

Experimental values for the isotopic cross-sections with their errors are listed in appendix A. The 54 Mn that 
could not be measured was obtained by a smooth interpolation between the neighboring isotopes so the value 
given in the tables is followed by (Interp.). This value is used to obtain integrated quantities as the mass or 
charge distributions and in the evaluation of the total reaction cross section also given in appendix A. 



Final results of the mean recoil velocity and the width of the velocity distributions for the various residual 
nuclei are presented in appendix B. Errors quoted here are due to the velocity reconstruction procedure above 
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described and to the magnetic-rigidity determination. In the case of a truncated velocity distribution, results 
partially interpolated are followed by (I) . The minus sign means that the recoil velocity is opposite to the original 
direction of the iron beam or in other words in the direction of the proton motion in the iron at rest system. 
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FIG. 4: Transmission factor as a function of the mass number of the residue for the five energies presented in this work 
(see text). 



III. RESULTS 



A. Isotope production cross sections 



Using the experimental method described above it was possible to measure at five different energies most of 
the residues produced in the spallation reaction of iron with cross-sections larger than 10 -2 mb, from cobalt 
(Z=27) down to oxygen (Z=8) or neon (Z=10) depending on the energy. At 1 A GeV, cobalt isotopes have not 
been measured. 

Figures [H O [3 [H and [5] show the isotopic distribution cross sections at the five beam energies. Error bars do 
not appear as they are smaller than the data points. The position of the maximum of these isotopic curves is 
correlated with the excitation energy transferred in the collision between the projectile and the target. In the 
case of a peripheral collision, in which the excitation energy is limited, only a few particles are evaporated by 
the fragment, leading to the population of isotopes close to stability. For more central collisions, the deposited 
excitation energy is larger and more neutron-deficient isotopes are produced due to the evaporation phase which 
favors the emission of neutrons. However, the tendency towards neutron-deficient isotopes is weaker than what 
is generally observed in heavy systems since, for iron, the Coulomb barrier is much smaller and the neutron to 
proton ratio in the projectile is also smaller. 

Isotopic cross-sections can be summed to obtain mass or charge distributions. Figure [TU] presents the mass 
distribution of the spallation residues for the five energies of the iron beam analyzed in this experiment. The 
residues are produced with different cross-sections depending on the energy of the projectile. The general trend 
of the data is globally as expected. As the beam energy increases, the deposited excitation energy becomes more 
and more important, leading in average to a stronger evaporation of nucleons, and finally to lighter evaporation 
residues. This is reflected by the substantial rise of the light fragment cross-sections between 300 and 1500 MeV 
per nucleon. As the total reaction cross section is overall rather constant over the studied energy range, this is 
compensated by a decrease of the production cross sections of the heaviest evaporation residues with increasing 
energy. It appears that masses around 46-47 are produced with a cross section almost independent of the beam 
energy. 
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FIG. 5: Isotopic production cross-sections of fragments from the reaction 56 Fe + p at 1.5 A GeV as a function of mass 
number 
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FIG. 6: Isotopic production cross-sections of fragments from the reaction 56 Fe + p at 1.0 A GeV as a function of mass 
number 
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FIG. 7: Isotopic production cross-sections of fragments from the reaction 56 Fe + p at 0.75 A GeV as a function of mass 
number 
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FIG. 8: Isotopic production cross-sections of fragments from the reaction 56 Fe + p at 0.5 A GeV as a function of mass 
number 
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FIG. 9: Isotopic production cross-sections of fragments from the reaction 56 Fe + p at 0.3 A GeV as a function of mass 
number 
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FIG. 10: Mass distribution of the residual nuclei in the spallation reaction Fe + p at the five different beam energies. 
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B. Comparison with other experimental data 

1. Reverse kinematics 

The present data can be compared with the ones obtained by W. R. Webber and collaborators using the 
reverse kinematics method. Measurements were performed on either a thick CH2 target (from ~ 2glcm 2 to 
~ 6g/cm 2 ) subtracting the carbon contribution 0, 0, 0], or a liquid-hydrogen target (1.52g/cm 2 ) [!, 0| at 
SATURNE. In both cases, the fragments were detected with a telescope of scintillators and Cerenkov counters. 

The charge distributions of the spallation residues for several iron beam energies from 330 to 1615 A MeV 
0,H)0| have been measured down to Z around 15. In Fig.QTJ these results (histograms), at beam energies close 
to ours, namely 1512, 1086, 724, 520 and 330 A MeV are compared with the present cross-sections (symbols) 
summed over masses to obtain the charge distribution. The overall agreement is satisfying in terms of variation 
with energy and charge of the residue. A systematic dependence of the element cross sections with the parity 
of Z is consistently observed in both experiments. The deviation factor, i.e. the average ratio between the two 
experiments has been calculated and is shown in Table [IJ The cross-section for Z=24 at 1512 A MeV, for which 
the Webber value is much larger than the neighboring cross sections and inconsistent with a general trend, 
is excluded. At the three highest energies, it is perfectly compatible with the precisions of both experiments 
(5% to 20% for Webber et al. and 9% to 15% here). At 300 A MeV (330 A MeV), the discrepancy is larger 
but still acceptable considering the different energies (10%) of the two measurements. The highest value (1-28) 
for the deviation factor is found at 500 A MeV (520 A MeV). Although this could be caused by a particular 
experimental problem at this energy, it is still compatible within the respective errors, especially if one bears 
in mind that at low energy both errors are larger: in our case because of the large transmission correction and 
in the case of Webber because of corrections for secondary reactions. The same reasons could explain the fact 
that, for a given energy, the disagreement is increasing with decreasing Z values, as it can be seen in Fig. 1111 
Another argument is that if we plot charge-changing cross sections as a function of the beam energy for various 
charges, our results at 500 A MeV are ~10% below a smooth interpolation based on the other measured energies 
whereas the Webber values are ~20% above the interpolation. 



Energy/ A (MeV) 


300 


500 


750 


1000 


1500 


Deviation factor 


1.23 


1.28 


1.01 


0.89 


0.88 



TABLE I: Average ratio of the charge- changing cross-sections measured by Webber et al. 0, 0, 0/ divided by the values 
from this experiment. 



The isotopic production cross-sections have also been measured previously but only at one energy (573 
A MeV), using a liquid-hydrogen target 0,0] and were limited to rather large cross-sections. The ratio between 
these values and the present data is displayed in Fig. 1121 including the respective errors. The lines represent 
the ratios of the cross sections at 573 A MeV and 500 A MeV computed with the INCL4-GEMINI combination 
of models. It shows that the difference in energy between the two experiments is not negligible for the lightest 
fragments, for which it can lead to differences of 30 to 40%. The agreement between the data is quite good for 
residues close in mass to iron but the difference increases for lighter isotopes. The value of the ratio is frequently 
hardly compatible with the expected value given by the line. 

Actually, one would expect a smooth variation of the mean value and of the width of the isotopic distribution 
with element charge. In Fig. [T3] are represented the mean mass-over-charge ratio as the function of Z, summing 
only the isotopes measured by both experiments. Clearly these quantities are more fluctuating in the Webber et 
al. experiment, in particular for potassium (Z=19) data and to a smaller extent for argon (Z=18) and titanium 
(Z=22) ones. The use of our full isotopic distributions, which extend much beyond the ones of Webber et al., 
does not make a large difference. 



2. Direct kinematics 

Results in direct kinematics have been obtained by R. Michel and collaborators 0, [34|, [3f| by irradiation of 
natural iron targets at different proton beam energies, allowing the determination of excitation functions from 
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FIG. 11: Nuclear-charge distribution of the residual nuclei for the five energies with scaling factors (2/1/0.5/0.25/0.125 
respectively from 1500 A MeV to 300 A MeV) applied for clarity. Points correspond to the present data, and solid 
histograms are data from Webber et al. @,Q[3] at close energies : 1512, 1086, 724, 520 and 330 A MeV. 
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FIG. 12: Ratio between isotopic cross-sections measured by Webber et al. at 573 A MeV 3, 4] and the present experiment 
at 500 A MeV for each element as a function of the mass number. Lines are theoretical predictions from INCL4-GEMINI 
for this ratio. 



a few tens of MeV to about 2 GeV. Some of the produced residual nuclei have been measured and identified by 
their gamma-ray decay spectrum or by mass-spectrometry. These data are compared to our experimental data 
in Fig. [M] and Fig. [15] Results can be split into "cumulative" and "independent" nuclei meaning that they are 
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z 

FIG. 13: Ratio of the mean nuclear-mass value to the charge for each element measured by Webber et al. at 573 A MeV 0, 
[5[ (black circles) and in the present experiment at 500 A MeV (open squares). Black squares are for the same quantity 
evaluated from all isotopes measured in this experiment. 

or not populated by a decay chain. In the case of cumulative cross-sections, our own cross-sections have been 
summed along the decay chain before comparing with Michel's data. The following isotopes: 36 Cl 42 K, 46 Sc, 
48 Sc 54 Mn and 52 Fe are "independent". 

Our data follow quite well, in most of the cases, the dependence on energy obtained in the Michel et al. 
experiment. This is very satisfying if we consider the difference between the two experimental methods. Some 
of the important differences that can be noticed may be due to the use of natural iron in the case of Michel's data. 
For instance, the observed higher cross-section for 52 Fe could come from a contribution of the (p,2n) reaction 
on 54 Fe adding to the (p,4n) on 56 Fe. Although there is only 6% of 54 Fe in natural iron, the effect should 
be non-negligible since (p,2n) is 40 times more probable than (p,4n) as deduced from our results. Conversely, 
the lower cross-sections found by Michel for 52 Mn and the higher one for 56 Co at high energy do not seem 
compatible with the tendency deduced from our isotopic distributions. 
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FIG. 14: Excitation functions of some residual nuclei produced in the spallation reaction of proton on iron. Open dots 
are the data of R. Michel et al. obtained by a direct irradiation 34. 35] and black triangles correspond to the present 
experimental data at 5 energies. Independent isotopes are indicated (Ind.). 
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FIG. 15: Excitation functions of some residual nuclei produced in the spallation reaction of proton on iron. Open dots 
are the data of R. Michel et al. obtained by a direct irradiation 34. 35] and black triangles correspond to the present 
experimental data at 5 energies. Independent isotopes are indicated (Ind.). 

Finally, we can say that the present results are qualitatively in good agreement with former measurements. 
The fact that we have complete isotopic distributions extending down to lighter nuclei than previously measured, 
on a wide range of energy, allows us to check the consistency of our own results and detect possible inconsistencies 
in other sets of data. 



C. Comparison with parametric formulas 

Since 1950, parametric formulas have been developed by astrophysicists with the aim of predicting the pro- 
duction cross-sections of the residual nuclei in spallation reactions. These formulas are used in case of light 
and intermediate nuclei present in the composition of the cosmic-rays like iron. In this section we present the 
comparison of our new experimental data with the results of three of these parametric formulas: Webber (3p| . 
EPAX 37] and Silberberg and Tsao [HH formulas. 
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1. Webber's formula 

This parametric formula has been developed by Webber et al. jj36| from the experimental data shown in the 
previous section. It is used in case of light spallation residues with Zi < 28 and for energies of the projectile 
E > 200 MeV. 

The form of this formula is : 



a(Ai, Z h E) = a {Zi, Z t ) ■ frfa, A h Z u A t ) ■ f 2 {E, Z h Z t ) 
for residual nuclei {Zi, Ai) of the spallation reaction on a target nuclei {Z t ,A t ) at energy E. 

• The first factor ao{Zi, Z t ) gives the charge distribution of the residues 

• f\{Zi, Ai,Zt, At) describes the isotopic curves (from their data at 573 MeV per nucleon) 

• f L 2{E, Zi,Z t ) gives the energy dependence 



Fe + p - Webber 




FIG. 16: Comparison between present results at 5 energies (symbols) and the results obtained with Webber's formula 
(solid lines). Scaling factors (2/1/0.5/0.25/0.125 respectively from 1500 A MeV to 300 A MeV) are applied for clarity. 



In Fig. [16] the comparison of our mass distribution with the predictions of the Webber's formula is shown. 
A rather good agreement is obtained at all the energies considered here for the heaviest residues, which are 
precisely those already measured by Webber et al. and used to determine the parameters of the formula. 
However, there seems to be some oscillations in the cross-sections not observed in the data. Actually, the 
charge distribution, not shown, is more accurately predicted by the formula. This comes from the fact that the 
isotopic distributions predicted by the Webber's formula have smaller widths than those obtained experimentally 
(see Fig. 1191) . A probable explanation is that only very few isotopic data were existing at the time when the 
formula was established. Therefore, the isotopic dependence could not be properly determined. Furthermore, 
the extrapolation of the parametric formula for light residues that are measured here for the first time shows 
an important discrepancy with the data. Even if this parametric formula can be useful for determining the 
production of the most produced spallation residues, this illustrates the danger of using parametric formulas 
outside the range in which they were adjusted. 
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2. The EPAX formula 

Epax formula [37l | was created with the aim of describing the production of residues in fragmentation reactions 
between heavy ions in what is call the limiting-fragmentation regime in which the residue production cross- 
section does not depend anymore on the energy of the projectile. Although not fully valid for protons at 
these energies, it might be instructive to know how close its predictions are to the present data. The limiting- 
fragmentation regime for the spallation reaction Fe + p is expected to be reached for energies of a few GeV per 
nucleon so here we can just expect the 1.5 A GeV data to be comparable with it. 

It can be used for spallation reactions with protons in the case of target nuclei of masses 18 < At < 187, 
although developed mainly for heavy-ion reaction . It is composed by two factors : 

a(Zi, Ai) = Ya ■ a(Z pro b - Z t ) 

with : 

• Ya a factor to describe the mass distribution of the fragments (Z i: A4) 



Y A = S 2 (A] /3 + AlH + Si) ■ P ■ exp[-P(A t - A,)} 

and InP — Pi- At +P1. Si, S2, Pi and P2 being adjusted parameters and A pro the mass of the projectile 
(one here for protons). 

• o~{Z pro b — Zi) describes the isobaric curves with Z pro t, as the charge for the maximal production. The 
various Z pro b values as a function of A defines the so called residue corridor in this approach. 

In Fig. [TTIour experimental results (symbols) are compared with the predictions of the EPAX formula. The 
experimental data at 1.5 GeV per nucleon are expected to be the ones closest to the limiting fragmentation 
regime, therefore we have renormalized the factor S2 so that the formula gives the total reaction cross-section 
measured at 1.5 GeV per nucleon (794 mb). Since the EPAX total cross section was 617 mb, this led to a 
multiplication by 1.28. 

It can be seen, as expected, that as the energy increases the mass distribution gets closer and closer to the 
EPAX prediction, with a quite good agreement at 1.5 GeV per nucleon. However, the lightest residues are still 
overestimated by the formula. The EPAX formula predicts also a more important evaporation of neutrons than 
seen in the isotopic cross-section data. In fact the measured N/Z ratio of the residues is higher than the one of 
the residue corridor which is used in the formula. 



3. Silberberg and Tsao 's formula 



The first version of this parametric formula has been developed in 1973 [40j with the experimental data 
measured by Rudstam [4l[ concerning the spallation residues in the spallation reaction p + Fe at 340 MeV. 
Various improvements, especially the beam-energy dependence, have been added in successive versions (H, l39j. 
It can be written as : 



a(A, Z, E)=a Q - f(A) ■ f{E) ■ e - p ^ AA ■ e -a\z-SA-TA>\> > Q . ^ . $ 

where : 

• (To is a normalization to the total reaction cross-section 

• f(A) and f(E) are factors used only in the case of target nuclei Z t > 30 

• e - p ( E ) AA represents the reduction in the production cross-section with the mass difference (A^4) between 
the residue and the target nuclei and an energy dependence through the P parameter 
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FIG. 17: Comparison between the present measured mass distributions (symbols) and the results obtained with the 
EPAX formula (solid line) . 

• e -R\ z ~SA-TA | describes the width and the position of the maximum in the isotopic and isobaric pro- 
duction 

• f2 takes into account the level structure of the residual nuclei 

• r\ is a factor for the pairing of protons and neutrons 

• $ represents an increase in the production of very light residues 

In Fig. [18] a comparison of this formula with the experimental results presented in this work is shown. In 
general, the agreement is very good for all energies between 10 % and 30 % at 300 MeV per nucleon where the 
discrepancy is larger. 

This last parametric formula appears as the most suitable to reproduce the present data, probably because 
of the largest data base used to derive it, which contains systems rather close to the ones studied here. These 
formulas are quite useful for quickly calculating production rates. Although some physical ingredients are 
present to derive them, more sophisticated approaches are needed to better handle the physics included in 
spallation reactions and to describe more fully other observables than cross-sections. 



4- Isotopic distribution shapes 



In the preceding sections only mass distributions were compared to the predictions of the parametric formula. 
It is also interesting to know how they reproduce the isotopic distributions. A powerful way to look at this is 
to compare the shape of the mass distributions of each element through the mean value and width of the mass- 
over-charge distributions as a function of Z. This is what is shown in Fig. fTS] in which the experimental results 
at 1500 MeV per nucleon (for better chance of agreement of EPAX) are compared with the three parametric 
formula. It can be seen that, as concerns the mean mass- over- charge, Webber's formula and EPAX agree rather 
well with the data while Silberberg-Tsao's predicts a slightly too high value. Regarding the widths, EPAX is 
acceptable and Webber tends to produce a too narrow mass distribution, maybe because the formula was fitted 
on his isotopic data which have a rather limited extension. Silberberg-Tsao gives a nearly constant width with 
Z, in contradiction to the experimental shape. This means that this formula, which gave the best agreement 
for mass distribution, should be used with caution if one wants to estimate isotope production cross-sections. 
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Fe + p - Silberberg and Tsao 




FIG. 18: Comparison between present experimental results at the five energies (symbols) and the results obtained with 
the Silberberg and Tsao's formula (solid lines). Scaling factors (2/1/0.5/0.25/0.125 respectively from 1500 A MeV to 
300 A MeV) are applied for clarity. 




FIG. 19: Mean values and width of the mass-over-charge distributions as a function of element charge measured at 
1500 MeV per nucleon compared with the predictions of the parametric formulas of Webber (solid line), Silberberg-Tsao 
(dashed-dotted line) and EPAX (dashed line). 



D. Comparison with models 

The design and optimization of spallation sources requires the knowledge of a large number of quantities 
directly related to spallation reactions in different materials and at various energies. Since exhaustive measure- 
ments of such a large amount of data are beyond the experimental possibilities, one needs to develop spallation 
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models with a good predictability that can be used in transport codes for simulations. This implies a deeper 
knowledge of the physics of the spallation reactions. 

Spallation is generally described by a two-step mechanism. The first stage of intra-nuclear cascade process 
(INC) governed by nucleon-nucleon collisions, leads to an excited nucleus after the ejection of a few energetic 
particles (p, n, it, d, a etc.). The second longest phase follows corresponding to the evaporative decay of the 
excited remnant nucleus with a possible competition with fission and Fermi break-up in some cases. Some 
approaches include also an intermediate stage of pre-cquilibrium to account smoothly for the transition to the 
full thcrmalization of the evaporating nucleus. 

Old INC models are still currently used in the high-energy transport codes employed for applications as 
Bertini [12] or ISABEL (43J) models. However, recently, a renewed interest for INC codes has been triggered 
by new available spallation data. Among them one could cite recent improvements on the INC codes found 
in phfl , |45| . In the present work, we have compared the experimental results of the spallation residues on iron to 
the predictions of three INC codes: Bertini code, ISABEL and INCL4 [44| . Since a long time, the first two ones 
are available in transport codes like LAHET3 (4|| and MCNPX [13] for simulations of macro-systems, Bertini 
(with pre-equilibrium) being used by default. INCL4 was only recently implemented in these code systems as 
well as CEM 45]. The basic physical assumptions are rather similar but differ in their implementation, for 
instance the way to develop the NN series of interactions, the way to treat Pauli blocking or the criterium to 
stop the INC stage. Note that we have used the implementation of ISABEL in LAHET3 which is blocked above 
1 GeV. But this does not means that this cascade is not valid at higher energies. 

For the second stage of the reaction, the most commonly used de-excitation code (and default option) in 
LAHET and MCNPX is the Dresner eva por ation code [HI complemented with the Atchison model for fission [3§] . 
It uses the Weisskopf-Ewing formalism [50] for the treatment of the evaporation, as do the more recent models 
ABLA [jyj and GEM (52j. Mainly, these three codes differ in the formulas and parameters used to described 
the level densities, the Coulomb barriers and the inverse reaction cross-sections. The Dresner code includes 
only the evaporation of light particles: neutrons, hydrogen and helium isotopes. The ABLA code has been 
mainly tuned for heavy systems with a particular interest on the fission description. In the version used in 
this work only neutrons, protons and alpha particles are evaporated. Furthermore shell and pairing effects as 
well as gamma decay were not taken here into account. The GEM code is a recent update of the Dresner 
model with new parameters and extends the Weisskopf-Ewing formalism to the evaporation of intermediate- 
mass fragments up to Z = 12. Actually the three codes (Dresner, ABLA and GEM) do not take into account 
in the evaporation process the angular momenta, which in fact are relatively small in spallation reactions 
induced by incident protons. Fission of heavy systems is described in a Bohr and Wheeler approach using 
phenomenological fragment distributions in Atchison and GEM. In ABLA fission is treated as a dynamical 
process taking into account the nuclear viscosity, and the fragment distribution is essentially obtained through 
the calculated population of states above the mass-asymmetric conditional saddle point. 

As will be shown in the following, conventional Weisskopf-Ewing evaporation may be not sufficient to account 
for our data. This is why we will also compare our results to models predicting the emission of intermediate- 
mass fragments through other mechanisms. The GEMINI model [H| treats evaporation of light particles within 
the Hauser-Feshbach formalism [54j . taking explicitly into account the angular momentum. Following the idea 
of Moretto (5f| that there should be a continuous transition between evaporation and fission, for all systems 
including light ones, the emission of intermediate fragments is handled as asymmetric fission in the Transition 
State Model. The transition between Hauser-Feshbach evaporation and asymetric fission can be chosen through 
a parameter: in the present work, this parameter has been set so that the Transition State Model is used for 
fragments above helium. Several other options exist in the code. We have used the ones recommended by the 
author. Some tries to vary them, although not exhaustive, do not reveal strong differences in the description of 
the present data. 

The SMM code is a numerical implementation of the Statistical Multifragmentation Model from [5(| often 
used to describe heavy-ion collisions in which multifragmentation is more li kely to arise. The parameters to 
describe the multifragmentation process are the standard ones as described in [571 ]. In particular, the asymptotic 
freeze-out volume is three times the initial one. The evaporation is treated in the Weisskopf-Ewing formalism 
up to fragment mass 18, and the lightest primary fragments decays arc treated by the Fermi break-up [58]. 

In the comparison between experimental data and model predictions, it is always difficult to disentangle the 
respective roles of the intra-nuclear cascade, which determines the characteristics of the remnant nucleus (charge, 
mass, angular momentum and excitation energy) at the end of the cascade stage, and of the de-excitation model. 
For instance, the under-prediction by the INCL4-ABLA combination of models of the light evaporation residue 
cross-sections observed for heavy systems [TH, H3] could be ascribed either to a too low excitation energy given 
by INCL4 or to a deficiency of ABLA at the highest excitation energies. However, some observables can be 
found that are more sensitive to one reaction stage or the other. In the following, we will try, as far as possible, 
to disentangle the influences of the intra-nuclear cascade and of the de-excitation stage in the comparison with 
the different observables. 
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1. Total reaction cross-section 



The total reaction cross-section is clearly one of the observables that depends only on the INC model since 
it is mainly related to the probability that the incident nucleon makes a collision with one nucleon of the 
target and that this collision is not blocked by the Pauli principle. In Fig. [201 we present the total reaction 
cross-sections obtained for the five energies analyzed in this work. They were calculated by summing up the 
isotope productions tabulated in the appendix. The summation has been done down to Z=8-10 depending on 
the bombarding energies. The contribution of the unmeasured isotopes have been estimated to be at most a few 
percents, i.e. smaller than the error bars. The fact that the lightest fragments could come from binary breakups 
and therefore leads to a possible double counting in the total reaction cross-section is also negligible. Actually, 
the two contributions play in opposit directions and even more or less compensate. Previous experimental data 
from the Barashenkhov compilation (59l | are also shown on this figure. A reasonable agreement is observed 
between most of the previous data and the present ones for both the absolute values and the behavior with 
the incident energy. The predictions of all the three INC codes agree with the data within the experimental 
accuracy, the difference between them being at most 10 %. This is not surprising since these INC models are 
known to generally well reproduce the total reaction cross-sections at energy above a hundred MeV 0, |46| . 
This observable cannot be used to discriminate between these three codes. 
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FIG. 20: Total reaction cross-sections of protons on iron as a function of the bombarding energy. Our five experimental 
data are compared to the compilation of previous experimental data from Barashenkhov [5£| and the values given by 
the three INC codes : Bertini, ISABEL (not available for E > 1 GeV in LAHET3) and INCL4. 



2. Mass and charge distributions 



In this section we examine the various model predictions compared to the mass or charge distributions 
obtained by summing the measured isotopic cross-sections. For completeness the light fragment cross sections 
analyzed in [2{| and obtained during the same experiment are also included at 1 GeV per nucleon. 

We first investigate the influence of the choice of the INC model. In Fig. [21] the mass and charge distributions 
of the residual nuclei produced at 1 GeV is shown and compared to the Bertini intra-nuclear cascade (plus 
pre-equilibrium) followed by the Dresner evaporation. Both mass and charge distributions lead to the same 
conclusions. The production yields of residues close to iron which are the major part of the spallation cross- 
section are underestimated while the yield of intermediate-mass residues is on the contrary overpredicted. The 
same conclusions were already obtained for heavy nuclei 18]. This behavior could be ascribed to a too high 
excitation energies at the end of the Bertini intranuclear cascade even after the introduction of a preequilibrium 
phase. A comparison is also shown with INCL4 followed by the same evaporation code. The calculations now 
predicts less excited remnants and a more satisfactory agreement is obtained for the heaviest residues but the 
light ones are still underestimated. It can be also noticed that, in both cases, the production of very light 
fragments is by far underpredicted. Another comparison is shown in Fig. [2U between the mass distribution 
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FIG. 21: Mass distribution of the spallation residues of iron at 1 A GeV compared to the Bertini and INCL4 INC models 
combined with the Dresner evaporation code. 



of the spallation residues and the predictions of two different INC models, now ISABEL and INCL4, followed 
by the ABLA evaporation. This last combination has been shown to reproduce satisfactorily many spallation 
data [44] in a wide domain and without adjustment of parameters. Both calculations give a similar good 
descriptions of the residues close to iron and an underprediction of intermediate and light nuclei cross-sections. 
This underprediction of INCL4- ABLA is in fact consistent with light evaporation residue cross-sections obtained 
from heavier nuclei (lead and gold) [44[ • Actually, for the heaviest nuclei which are mainly formed in peripheral 
collisions with low excitation energy, the evaporation plays a less important role than the intra-nuclear cascade 
since only a very little number of nucleons is evaporated. The fact that both calculations have the same behavior 
and are rather good for heavy residues suggests that the underprediction of the light residues is not due to a 
lack of excitation energy. Indeed, we have seen in the comparison with Bertini in Fig.[5T]that a larger excitation 
energy does lead to a larger production of light fragments but to the detriment of heavy ones which cannot be 
counterbalanced by playing with evaporation models. This rather indicates that the problem comes from the 
de-excitation stage. In the following we will not consider anymore the Bertini model for which many shotcomings 
have been pointed out [l3| , [l8[ , [6(| ■ We will mainly restrict the comparisons with various de-excitation models 
using INCL4 in the first stage, since ISABEL generally gives similar results. 

Figure |2"31 shows the INCL4 intranuclear cascade coupled with the GEM model, which takes into account also 
evaporation of intermediate-mass fragments. The calculated cross-sections for the intermediate mass residues 
are improved comparatively to ABLA. However one observes a slight underestimation of the residues close to 
iron and the underprediction of the very light fragments still persists for masses slightly smaller than with 
ABLA. 

From the comparison with the three evaporation models (Dresner, ABLA and GEM) and the remark con- 
cerning excitation energy made above, it can be presumed that standard evaporation models, even including 
the emission of IMF (GEM), cannot reproduce the bulk of our data. This is why we tried other models which 
include other de-excitation modes. 

On Fig. [23l are also shown the predictions of GEMINI. If on the heavy fragments the results are slightly 
less satisfactory than with GEM, the behaviour for A lower than 30 is significantly improved. Probably due 
to its capability of predicting asymmetric fission in the Transition State Model prescription, GEMINI appears 
as the best suited code to reproduce the bulk of the data except at the lower energy (300 MeV per nucleon). 
Actually at 300 MeV per nucleon, all the calculations, whatever the choice of INC or de-excitation models, start 
to deviate from experiment around A equal 48. 

Even-odd disymmetry of the cross sections are clearly visible on an unlarged picture of the Z distribution 
at 1 GeV per nucleon in Fig. [M] representative also of other energies. In spite of a small underprediction of 
the absolute cross-sections with GEMINI, the ratios between odd and even Z cross-sections are very close to 
the experimental ones. Whereas GEM gives a too strong effect, ABLA predicts (with the present version) a 
slightly too small even-odd effect. But for the largest cross-sections above 18 the INCL4-ABLA remains the 
more precise prediction of the experimental values. 

Another mechanism that could be invoked to explain our large yields of light fragments could be the onset of 
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FIG. 22: Mass distribution of the spallation residues of iron at 1 A GeV compared to two different INC codes (INCL4 [44} 
and ISABEL [11]) combined with the ABLA evaporation code [Hit ]. 




A 

FIG. 23: Spallation residue cross-sections of iron as a function of there mass number compared with a calculation 
with INCL4 coupled with GEM (dashed lines) or GEMINI (continuous lines). Points are data of the present paper 
complemented for low masses at 1 GeV by the ones of [29| ] obtained during the same experiment. 

multifragmentation at the highest excitation energies [2^| . The coupling of INCL4 with the multifragmentation 
model, SMM, is shown in Fig. [55] The model well describes the heavy residues and the ones with mass between 
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FIG. 24: Charge distribution of the spallation residues of iron at 1 A GeV compared to the INCL4 cascade coupled to the 
deexcitation models GEM, GEMINI and ABLA. 



20 and 30. However, it overpredicts the lighest fragments and disagrees strongly with the data in the A region 
30-45. The contribution of fragments produced by multifragmentation is shown as the dashed curve in the figure 
(multifragmentation events being identified by the entry into the multifragmentation routine in the code |6l|). 
The major part of the light fragment cross-section is produced by multifragmentation while masses above 25 
are mostly originating from evaporation. However, it is likely the opening of multifragmentation that causes 
the hole in the region A = 30 — 45, not observed experimentally. Our results are at variance with what was 
found in [29], where SMM coupled to another INC model (from [62]) was giving a good agreement with the 
data, provided that a pre-equilibrium stage was added. With INCL4, which as explained in j44j, handles what is 
often called the pre-equilibrium stage, the best agreement with the whole set of data is obtained with GEMINI. 




FIG. 25: Calculation of the cross sections with INCL4 and SMM (solid line) compared to data points at 1 A GeV. The 
multifragmentation contribution to the calculation is the dashed line. 



However, a clear conclusion on the mechanism responsible for the light and intermediate fragment production 
is difficult and would need more constraining information. It seems rather clear that the traditional Weisskopf- 
Ewing evaporation as used in ABLA or even in GEM, which evaporates IMFs, miss the production of the lighest 
nuclei. However, the reason for the success of GEMINI, Hauser-Fesbach treatment or asymmetric fission from 
the Transition-state-model, is not fully understood and a possible contribution of multifragmentation is not 
ruled out. Forthcoming exclusive experiments will probably help to clarify the situation by an identification of 
the various fragments emitted in coincidence during the de-excitation stage of the reaction. 
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3. Isotopic distributions 

In this experiment, more than 500 individual isotopic cross sections have been measured which have been 
compared systematically to calculations done with the four different de-excitation codes (ABLA, GEM, GEMINI 
and SMM) coupled with INCL4. As an example, the comparison of GEMINI (full line) and ABLA (dashed 
line) with a selection of measured isotopic cross sections at 1 A GeV is shown on Fig. [23 Except for the better 
level of cross sections for light residues from GEMINI, already seen when looking at the mass distributions, it 
is difficult to conclude about the detailed quality of each model. 




20 40 60 20 40 60 20 40 



A 

FIG. 26: A selection of isotopic distributions of cross-sections measured at 1 A GeV compared with INCL4-GEMINI 
calculations (continuous lines) and INCL4-ABLA (dashed lines). 

A more powerfull way to make the comparison is to look at the shapes of the isotopic distributions for each 
element through the mean atomic mass ({A}) divided by the charge Z of the element and the width (root mean 
square) of the measured (or computed) distributions. Figure [2"T1 presents the comparison of these quantities with 
GEMINI at the five energies while Fig. [28] shows the results INCL4 coupled to either ABLA, GEM or SMM 
at 1 A GeV. Actually, it is remarkable that the deviations between models and experiment are qualitatively 
independent of the beam energy. This can be checked on Fig. [57] for GEMINI but holds also for the comparison 
with the other models. For this reason, the comparison with the other three models is shown only at 1 A GeV 
in Fig. [28] But again, the following conclusions are the same at all the energies. 

For Z equal 25, 26 and 27 (not measured at 1 A GeV), cross sections are dominated by the cascade, leaving 
the remnant nucleus with very little excitation energy. Therefore, the choice of the evaporation model play 
practically no role and basically the (A) jZ is perfectly reproduced. The average value of the isotopic distribution 
{(A)/Z) is actually very well predicted by GEM and GEMINI on all the range (down to Z equal 8 to 9), with 
the correct odd-even effects, whereas ABLA gives a value systematically too small. The SMM model gives a 
correct centroid down to Z equal 20 but is the worst below this value with a distribution centered one mass 
below the data at lower Z. As regards the width of the distributions, none of the models is good on all the Z 
range. The widths computed from GEMINI are systematically a little too wide. With GEM and ABLA, they 
are too wide only in the range 20-25, otherwise they are very close to the data. For SMM, it is the contrary, 
rather good at high Z but too narrow for lower charges. This fact was already noticed in [631 ]. 

All this shows that none of the de-excitation models is perfect. However, taking into account the information 
from both the cross-sections and the isotopic distribution shapes, it can be concluded that the GEMINI gives 
the best agreement with our data. 

4- Recoil velocities. 

Concerning the kinetic characteristics of the fragments, we show in Fig. [29] a comparison between the ex- 
perimental mean longitudinal recoil velocities for each mass and the ones calculated with the INCL4 model 
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FIG. 27: Mean atomic mass A over Z ({A)/Z) and rms (a(A)) of the isotopic cross section distributions as a function 
of the atomic number Z, at the five bombarding energies, compared to calculations done with the intranuclear cascade 
INCL4 model coupled with GEMINI. The calculated values have been averaged over the actually measured isotopes. 



combined with ABLA or GEMINI at 1 A GeV. The same comparison is also done for the width (root mean 
square) of the longitudinal distribution (right part of the figure) . 

One can observe an important discrepancy between the experimental mean recoil velocities and values pre- 
dicted by the models. It is worthwhile to note that the experimental data decrease much more slowly with 
decreasing mass than the values predicted by the codes. Furthermore, they seem to saturate at a mass value 
of 35. This saturation is not seen with ABLA. Only GEMINI shows a clear tendency towards saturation below 
A=30. For the widths, on the contrary, the agreement with the experimental data is better, especially when 
using GEMINI. This behavior, presented here at 1 GeV, is very similar at the others energies analyzed in this 
experiment (Fig. [311)1 . The better agreement with GEMINI could be due to the existence of binary decays 
in de-excitation phase that reduces the mean longitudinal velocity of the final residual nuclei since the recoil 
momentum will originate from a heavier nucleus and will be split between two partners emitted in an arbitrary 
direction with respect to the beam one. In the same figure, are also shown the predictions from systematics of 
Morrisey [64| . which more or less give the correct slope for large mass but miss the saturation observed in the 
data. Actually, the two other de-excitation models GEM and SMM, not shown here, give results rather similar 
to ABLA: rather good for the widths but a slope too steep and an inability to describe the saturation of the 
mean values. 

The fact that the mean recoil velocities for the heaviest masses is not well predicted cannot be ascribed to the 
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FIG. 28: Same as fig[27]but only at 1 A GeV for comparison with calculations done with the intranuclear-cascade INCL4 
model coupled with ABLA (full line), GEM (dashed line), and SMM (dashed-dotted line). 

de-excitation models but should rather raise questions on the intranuclear cascade. This is why we have also 
performed a calculation using ISABEL coupled to ABLA, which is presented in Fig. [2H Obviously, ISABEL 
better reproduces both the mean values and the widths for masses larger than 50, indicating a possible deficiency 
of INCL4 in the recoil velocity determination. Actually, a simila r sy stematic deviation of INCL4 concerning the 
mean velocities has already been noticed for Pb+p at 1 GeV/A [18|. However, the general trend of the ISABEL 
calculation on the whole mass range leads to the same conclusion that it is unable to give the correct slope and 
saturation effect of the experimental data. 



IV. CONCLUSION 

The spallation residues produced in the bombardment of 56 Fe at 1.5 , 1.0, 0.75, 0.5 and 0.3 A GeV on a 
liquid-hydrogen target have been studied using the reverse kinematics technique and the Fragment Separator 
at GSI (Darmstadt). This technique has permitted the full identification in charge and mass of all isotopes 
produced with cross-sections larger than 10 -2 mb down to Z = 8. Their individual production cross-sections 
and recoil velocities at the five energies have been obtained. 
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FIG. 29: Mean (left) and Root Mean Square (right) values of the longitudinal recoil velocity distribution for spallation 
residues at 1 A GeV versus their atomic mass. Open circles are the experimental values. Down triangles are predictions 
from the INCL4-ABLA calculation, up empty triangles from the ISABEL-ABLA and up full triangles from the INCL4- 
GEMINI calculation. Velocities are expressed in the beam ( 56 Fe) rest frame and with a minus sign as being opposite to 
the iron beam direction. 



Fe+p, INCL42 + GEMINI 




Fe+p, INCL42 + GEMINI 



0.5 


0.5 


0.5 



^ 
>f 0.5 

S o 

Sh 

0.5 



1 




: 1500 A MeV 

.... 1 .... 1 ... . 


.... 


. . . . i T^T . 


| % 




: 1000 A Mt \ 

I I I , i , I , , i I I I I 


, , , , i , , , , i , , , , 


■ 1 °io D 




I °°o 

: 750 A MeV 
.... i .... i ... . 


.... 


. . . . i riZ . 


■ °Oc 




Cc 

: 500 A MeV 

, , , , i , , , , i , , , , 


, , , . 


. . . . i . y7 , 






! 0oc ° 

: 300 A MeV 
.... i .... i ... . 


°°°Oo 

.... 





10 20 



30 40 

A 



50 60 



FIG. 30: Mean values and r.m.s. of the longitudinal recoil velocity distribution for spallation residues versus their atomic 
mass at all beam energies. Open circles are experimental values. Down triangles are predictions from the INCL4 + 
GEMINI calculation. The lines are the Morrissey systematics . 



The production cross-sections have been compared with the previously existing data, either charge-changing 
cross-sections with a few isotopic cross-sections at one energy measured in reverse kinematics or excitation 
functions for a limited number of isotopes obtained by 7-spectrometry in direct kinematics. Globally, our 
results were found in good agreement with former data. This comparison also showed that our experimental 
method leads to a much more complete picture of the residue production than what was possible before with 
the few scattered results, allowing sometimes to detect possible inconsistencies in other sets of data. 

Comparisons with parametric formulas, often used in astrophysics, have been performed: the Webber formula 
gives rather good predictions of the charge distributions but produces too narrow isotopic distributions. It also 
totally fails for the lightest nuclei (below A equal 30-35) in the region not measured at the time when this 
formula was derived. The EPAX formula (once renormalized to give the correct total reaction cross-section) is 
usable only in the limiting fragmentation regime, apparently not yet fully reached at 1.5 A GeV. However, it 
nearly gives the right A dependence of the cross sections at our highest energy. The best formula seems to be 
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the Silberberg and Tsao one, which is in very good agreement with the experimental mass distributions and 
mean value of the isotopic distributions at all the energies except at 300 MeV (as all the models). The use 
of parametric formulas can be of great help for a fast estimation for certain applications, but the example of 
Webber's illustrates the possible danger of using parametric formulas outside the range on which they have 
been adjusted. Our data could certainly be used to derive new, more relaible parametric formulas for use in 
cosmic-ray propagation codes. 

Predictions of different intranuclear-cascade models (Bertini, ISABEL and INCL4) combined with different 
de-excitation models (Dresner, ABLA, GEM, SMM and GEMINI) have been confronted to the new experimental 
data. INCL4 or ISABEL combined with standard Weisskopf-Ewing evaporation models as ABLA or GEM give 
a good description of the residual production close in charge to iron but they underpredict systematically the 
light evaporation residues in the mass region 20-30. This fact, together with the saturation observed in the 
experimental longitudinal velocity at low masses, could be an indication that another de-excitation mechanism 
has to be considered. A de-excitation including a possible contribution from multifragmentation, as treated 
by SMM, improves significantly the predictions of light and intermediate mass fragments but at the detriment 
of residues in the region A — 30 — 45. SMM also misses the saturation of the recoil velocity and do not 
properly predict the isotopic distribution mean values and widths. The best overall agreement with the data 
is obtained with GEMINI. This GEMINI model gives a rather precise account of all cross sections measured 
here as a function of the beam energy. The recoil velocities, although not perfect, are the closest to the 
experimental values and the mean values and widths of the isotopic distributions are rather well reproduced. 
Other authors 0, EH have found that generally GEMINI reproduces very well the energy spectra of both 
light charged particles and intermediate mass fragments in a wide range of incident energies and target masses. 
Similar conclusions (best agreement with GEMINI) have been reached by [65[ using as intra-nuclear model 
the Cascade-Exciton Model, coupled with GEM, GEMINI and SMM and compared with these data taken 
from the C. Villagrasa-Canton PhD [66|. In [2i|, with another INC coupled with a preequilibrium stage, the 
deexcitation code SMM was found to give the best agreement with the 1 GeV data. It is obviously difficult to 
definitively conclude on the production mechanism of the intermediate and light mass fragments and probably 
only additional experimental information on correlations between residual nuclei and light particles could bring 
aswers to the questions addressed here. 

As regards to the potential interest of the present data for applications, we supply isotopic cross-sections that 
can be used to directly estimate the change in chemical composition that could occur in an ADS window made 
predominantly of iron and recoil velocities to calculate damages due to atomic displacements (DPA) [66| . 
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VI. APPENDIX A: CROSS SECTIONS. 



1500 A McV a R (mb) = 822 ± 73 

Z A a(mb) Z A a(mb) Z A a(mb) 

27 54 0.035+0.0021 23 50 16.272+2.383 20 44 8.0592+0.798 

27 55 0.245±0.0104 23 51 7.179±0.446 20 45 1.8149±0.174 

27 56 0.413+0.0187 23 52 1.274+0.073 20 46 0.4263+0.039 

26 51 0.016+0.0021 23 53 0.251+0.012 20 47 0.0497+0.004 

26 52 0.309+0.0187 23 54 0.004+0.001 20 48 0.0052+0.001 

26 53 2.865+0.1338 22 43 0.093+0.009 19 37 0.0683+0.008 

26 54 22.204+0.8962 22 44 1.309±0.129 19 38 1.2762+0.141 

26 55 58.838+2.2010 22 45 7.067+0.675 19 39 8.4339+0.918 

25 49 0.018+0.0021 22 46 19.052+1.767 19 40 12.7248+1.364 

25 50 0.329+0.0249 22 47 20.863+1.941 19 41 10.0161+1.076 

25 51 3.612+0.2240 22 48 15.246+1.957 19 42 4.2599+0.443 

25 52 15.089+0.8523 22 49 4.051+0.325 19 43 1.4303+0.147 

25 53 37.982+1.7678 22 50 1.078+0.081 19 44 0.2596+0.026 

25 54 42.812(Intcrp.) 22 51 0.103+0.006 19 45 0.0434+0.004 

25 55 32.885+1.2297 22 52 0.011+0.001 19 46 0.0041+0.001 

24 46 0.002+0.0004 21 41 0.027+0.004 18 35 0.1706+0.020 

24 47 0.054+0.0052 21 42 0.918+0.097 18 36 2.7062+0.304 

24 48 0.724+0.0622 21 43 6.840+0.702 18 37 8.8348+0.985 

24 49 5.043+0.4032 21 44 14.639+1.447 18 38 14.1045+1.551 

24 50 18.522+1.3723 21 45 17.225+1.687 18 39 8.6291+0.947 

24 51 27.886+1.7288 21 46 8.934+1.004 18 40 3.0793+0.331 

24 52 31.954+4.2184 21 47 2.989+0.267 18 41 0.7825+0.084 

24 53 10.052+0.4695 21 48 0.572+0.049 18 42 0.1602+0.017 

24 54 3.848+0.1565 21 49 0.093+0.007 18 43 0.0207+0.002 

24 55 0.054+0.0021 21 50 0.006+0.001 18 44 0.0021+0.001 

23 45 0.046+0.0041 20 39 0.105+0.011 17 33 0.0837+0.009 

23 46 0.693+0.0642 20 40 2.115+0.227 17 34 1.6685+0.192 

23 47 5.372+0.4797 20 41 9.204+0.982 17 35 8.7775+1.001 

23 48 15.447+1.3179 20 42 16.363+1.695 17 36 10.6805+1.203 



23 49 23.728+1.9001 20 43 14.268+1.500 17 37 6.9602+0.815 



1500 A MeV 



Z A 


a(mb) 


Z A cr(m6) 


Z 


A cr(mb) 


17 38 


1 872+0 2066 


14 29 8 920+1 055 


11 


25 1.0696+0.130 


17 39 


536±0 0589 


14 30 8.572+1.042 


11 


26 0.1719+0.021 


17 40 


098+0 0103 


14 31 2 549+0 314 


11 


27 0.0309+0.003 


17 41 


018+0 0021 


14 32 626+0 073 


10 


19 0.2531+0.032 


16 31 


027+0 0031 


14 33 072+0 008 

1 I ' } ' ) V ' ■ V7 1 _ _i_ V / . V/ V/ V.! 


10 20 3.2281+0.402 


16 32 


3 679+0 4306 


14 34 0.015+0.002 


10 


21 3.1828+0.446 


16 33 


9.282+1.0728 


13 25 0.246+0.030 


10 


22 3.0172+0.389 


16 34 


11.562+1.3826 


13 26 2 751+0 331 


10 


23 0.5833+0.071 


16 35 


4 865+0 5648 


13 27 7 469+0 886 


10 24 0.1378+0.016 


16 36 


1 547+0 1745 


13 28 4 550+0 545 


9 


17 0.2251+0.029 


16 37 


259+0 0289 


13 29 2 090+0 325 


9 


18 1.8555+0.234 


16 38 


047+0 0052 


13 30 366+0 043 


9 


19 2.0827+0.301 


16 39 


007+0 0010 


13 31 078+0 009 


9 


20 1.7157+0.249 


15 29 


171+0 0206 


12 23 0.498+0.061 


9 


21 0.5777+0.072 


15 30 


2 126+0 2508 


1 2 24 4 837+0 589 


9 


22 0.0874+0.010 


15 31 


8 475+0 9959 


12 25 5 600+0 676 


8 


16 6.2764+0.799 


15 32 


7.057+0.8514 


12 26 4.813+0.583 


8 


17 2.3020+0.364 


15 33 


3.707+0.4303 


12 27 1.185+0.142 


8 


18 1.6292+0.270 


15 34 


0.851+0.0980 


12 28 0.315+0.038 


8 


19 0.2845+0.035 


15 35 


0.193+0.0217 


12 29 0.033+0.004 





0.0000+0.000 


15 36 


0.022+0.0021 


11 21 0.226+0.028 





0.0000+0.000 


15 37 


0.002+0.0010 


11 22 2.244+0.276 





0.0000+0.000 


14 27 


0.453+0.0547 


11 23 4.967+0.608 





0.0000+0.000 


14 28 


6.317+0.7539 


11 24 2.594+0.319 





0.0000+0.000 



1000 A MeV o R {mb) = 811 ± 76 

Z A a(mb) Z A a(mb) Z A a(mb) 

26 51 0.016+0.0010 23 52 1.259+0.071 20 44 8.6501+0.896 

26 52 0.333+0.0187 23 53 0.246+0.011 20 45 1.8676±0.178 

26 53 3.319+0.1556 23 54 0.004+0.001 20 46 0.4159+0.038 

26 54 18.977+0.7686 22 42 0.002+0.001 20 47 0.0476+0.004 

26 55 47.982+1.7913 22 43 0.071+0.007 20 48 0.0041+0.001 

25 49 0.017+0.0010 22 44 1.347+0.134 19 37 0.0455+0.005 

25 50 0.373+0.0280 22 45 7.999+0.764 19 38 1.1211+0.123 

25 51 3.971+0.2457 22 46 22.716+2.106 19 39 8.3398+0.909 

25 52 15.927+0.9010 22 47 24.143+2.162 19 40 14.1478+1.517 

25 53 39.583+1.8373 22 48 17.468+2.361 19 41 10.5654+1.169 

25 54 46.679(Intcrp.) 22 49 4.245+0.340 19 42 4.4129+0.496 

25 55 35.410+1.3189 22 50 1.162+0.086 19 43 1.3951+0.143 

25 56 0.338+0.0135 22 51 0.100+0.006 19 44 0.2565+0.026 

24 46 0.002+0.0010 22 52 0.010+0.001 19 45 0.0403+0.004 

24 47 0.048+0.0041 21 41 0.019+0.002 18 34 0.0031+0.001 

24 48 0.789+0.0674 21 42 0.882+0.091 18 35 0.1013+0.011 

24 49 5.925+0.4737 21 43 7.427+0.762 18 36 2.2245+0.250 

24 50 21.256+1.5765 21 44 17.367+1.717 18 37 8.4183+0.939 

24 51 32.184+1.9900 21 45 19.221+1.837 18 38 14.8374+1.633 

24 52 34.246+8.0637 21 46 10.876+1.470 18 39 8.6601+0.992 

24 53 10.341+0.4799 21 47 3.198+0.286 18 40 3.1145+0.334 

24 54 3.558+0.1430 21 48 0.586+0.050 18 41 0.7143+0.076 

23 44 0.002+0.0010 21 49 0.089+0.007 18 42 0.1478+0.016 

23 45 0.037+0.0031 21 50 0.005+0.001 18 43 0.0186+0.002 

23 46 0.736+0.0684 20 38 0.003+0.001 17 32 0.0031+0.001 

23 47 6.104+0.5450 20 39 0.074+0.008 17 33 0.0486+0.005 

23 48 18.284+1.5613 20 40 2.012+0.215 17 34 1.2491+0.145 

23 49 27.086+2.1654 20 41 9.679+1.033 17 35 7.8911+0.901 

23 50 19.590+3.6884 20 42 18.973+1.966 17 36 10.3138+1.195 



23 51 7.109+0.4393 20 43 16.390+2.018 17 37 6.7298+0.807 



1000 A MeV 



Z A 


a(mb) 


Z 


A a(mb) 


17 38 


1 608+0 1 767 


13 27 5.202+0.625 


17 39 


464±0 0506 


13 28 3.173+0.384 


17 40 


080+0 0083 


13 


29 1.393+0.171 


16 30 


004+0 0010 


13 


30 0.229+0.027 


16 31 


138+0 0165 


13 


31 0.047+0.005 


16 32 


2 611+0 3056 


12 


22 0.010+0.001 


16 33 


7 803+0 9025 


12 


23 0.158+0.020 


16 34 


11.850+1.6129 


12 


24 2.482+0.309 


16 35 


4 211+0 4977 

I ■ — 1 1 i_ V ' . 1 - / 1 ( 


12 25 3.793+0.458 


16 36 


1.270+0.1425 


12 


26 3.160+0.381 


16 37 


0.202+0.0227 


12 


27 0.760+0.092 


16 38 


037+0 0041 


12 28 0.178+0.022 


15 28 


004+0 0010 


11 


20 0.009+0.001 


15 29 


078+0 0093 


11 


21 0.087+0.011 


15 30 


1 355+0 1 600 


11 


22 1.050+0.132 


15 31 


6 668+0 7843 


11 


23 3.091+0.389 


15 32 


6 284+0 9288 


11 


24 1.480+0.222 


15 33 


3 066+0 3622 


11 


25 0.629+0.076 


15 34 


623+0 0722 


11 


26 0.096+0.011 


15 35 


1 37+0 01 55 


10 


18 0.007+0.001 


14 26 


010+0 0010 


10 


19 0.140+0.019 


14 27 


0.180+0.0217 


10 


20 1.381+0.178 


14 28 


3.852+0.4600 


10 


21 2.199+0.280 


14 29 


6.525+0.7746 


10 22 1.664+0.210 


14 30 


7.040+0.8333 


10 


23 0.299+0.038 


14 31 


1.870+0.2269 


9 


17 0.072+0.009 


14 32 


0.431+0.0505 


9 


18 0.706+0.094 


14 33 


0.047+0.0052 


9 


19 1.452+0.196 


13 24 


0.006+0.0010 


9 


20 0.891+0.121 



13 25 0.094+0.0113 9 21 0.300+0.037 



13 26 1.449+0.1742 8 17 1.585+0.211 



750 A MeV a R (mb) = 767 ± 66 

Z A a(mb) Z A a(mb) Z A a(mb) 

27 54 0.035+0.0021 23 50 19.689+2.487 20 44 8.3324+0.831 

27 55 0.333+0.0145 23 51 7.647+0.473 20 45 1.8625±0.178 

27 56 0.474+0.0228 23 52 1.270+0.073 20 46 0.3901+0.036 

26 51 0.011+0.0010 23 53 0.225+0.010 20 47 0.0445+0.004 

26 52 0.300+0.0176 23 54 0.003+0.001 20 48 0.0041+0.001 

26 53 3.437+0.1628 22 43 0.045+0.004 19 37 0.0207+0.002 

26 54 20.841+0.8433 22 44 1.016+0.100 19 38 0.7136+0.079 

26 55 52.837+1.9738 22 45 7.983+0.763 19 39 7.2694+0.792 

25 49 0.012+0.0010 22 46 22.428+2.081 19 40 11.9543+1.283 

25 50 0.298+0.0228 22 47 25.267+2.274 19 41 9.2923+1.138 

25 51 4.366+0.2717 22 48 18.309+1.823 19 42 3.6528+0.384 

25 52 16.830+0.9529 22 49 4.577+0.367 19 43 1.3052+0.134 

25 53 42.696+1.9835 22 50 1.145+0.085 19 44 0.2306+0.023 

25 54 47.052(Intcrp.) 22 51 0.093+0.006 19 45 0.0372+0.003 

25 55 34.496+1.2888 22 52 0.009+0.001 19 46 0.0031+0.001 

25 56 0.429+0.0197 21 41 0.011+0.001 18 35 0.0413+0.005 

24 47 0.034+0.0031 21 42 0.619+0.064 18 36 1.3831+0.156 

24 48 0.653+0.0560 21 43 7.052+0.724 18 37 6.7799+0.757 

24 49 6.143+0.4923 21 44 16.154+1.598 18 38 12.3038+1.355 

24 50 22.252+1.6501 21 45 19.447+1.895 18 39 7.2843+0.817 

24 51 34.636+2.1424 21 46 9.974+1.004 18 40 2.4932+0.269 

24 52 35.666+6.1359 21 47 3.387+0.302 18 41 0.5592+0.060 

24 53 10.957+0.5099 21 48 0.575+0.049 18 42 0.1034+0.010 

24 54 3.408+0.1379 21 49 0.087+0.007 18 43 0.0155+0.002 

24 55 0.037+0.0021 21 50 0.005+0.001 17 33 0.0196+0.002 

23 45 0.023+0.0021 20 39 0.039+0.004 17 34 0.7211+0.084 

23 46 0.569+0.0528 20 40 1.394+0.149 17 35 5.9777+0.686 

23 47 6.337+0.5657 20 41 8.992+0.961 17 36 7.7320+0.873 

23 48 18.897+1.6131 20 42 16.725+1.741 17 37 5.0809+0.578 



23 49 30.515+2.4399 20 43 15.585+1.707 17 38 1.2656+0.141 



750 A MeV 

Z A a(mb) Z A a(mb) 

17 39 0.322±0.0351 13 28 1.841±0.225 
17 40 0.052±0.0052 13 29 0.714±0.087 
16 31 0.046±0.0062 13 30 0.088±0.010 
16 32 1.437±0.1683 13 31 0.020±0.002 
16 33 5.539±0.6412 12 23 0.052±0.008 
16 34 8.110±0.9386 12 24 1.193±0.147 
16 35 2.919±0.3366 12 25 1.968±0.245 
16 36 0.896±0.1012 12 26 2.149±0.277 
16 37 0.113±0.0124 12 27 0.538±0.071 
16 38 0.022±0.0021 12 28 0.071±0.008 
16 39 0.003±0.0010 12 29 0.002±0.001 
15 29 0.025±0.0031 11 21 0.030±0.005 
15 30 0.664±0.0795 11 22 0.526±0.085 
15 31 4.572±0.5449 11 23 1.810±0.251 
15 32 4.348±0.6440 11 24 0.779±0.100 
15 33 2.074±0.2477 11 25 0.249±0.031 
15 34 0.368±0.0433 11 26 0.058±0.007 
15 35 0.074±0.0083 10 20 0.638±0.103 
15 36 0.005±0.0010 10 21 1.216±0.183 
14 27 0.053±0.0072 10 22 0.968±0.166 
14 28 1.964±0.2341 10 23 0.047±0.007 
14 29 3.960±0.4765 10 24 0.014±0.002 
14 30 4.510±0.5415 9 18 0.342±0.057 
14 31 1.062±0.1289 9 19 0.897±0.162 
14 32 0.239±0.0278 9 20 0.469±0.088 
14 33 0.013±0.0021 9 21 0.097±0.014 
14 34 0.004±0.0010 8 17 0.880±0.149 
13 25 0.030±0.0041 8 18 0.431±0.083 
13 26 0.702±0.0866 0.000±0.000 



13 27 3.243±0.4247 0.000±0.000 



500 A MeV a R (mb) = 660 ± 53 

Z A a(mb) Z A a(mb) Z A a(mb) 

27 53 0.003+0.0010 23 49 27.507+2.593 20 43 11.4146+1.372 

27 54 0.059+0.0021 23 50 19.027+2.696 20 44 6.7186+0.745 

27 55 0.485+0.0197 23 51 6.266+0.387 20 45 1.2635+0.121 

27 56 0.611+0.0248 23 52 0.949+0.054 20 46 0.2558+0.024 

26 51 0.020+0.0010 23 53 0.152+0.007 20 47 0.0268+0.002 

26 52 0.416+0.0238 23 54 0.002+0.001 20 48 0.0021+0.001 

26 53 3.855+0.1789 22 43 0.052+0.005 19 37 0.0227+0.003 

26 54 23.507+0.9454 22 44 1.052+0.104 19 38 0.5010+0.056 

26 55 56.647+2.1110 22 45 6.650+0.635 19 39 4.4844+0.489 

25 49 0.019+0.0021 22 46 20.458+1.897 19 40 8.0008+0.859 

25 50 0.387+0.0289 22 47 21.553+2.204 19 41 6.3328+0.763 

25 51 4.490+0.2771 22 48 17.702+2.230 19 42 2.7133+0.313 

25 52 18.309+1.0328 22 49 3.529+0.282 19 43 0.7948+0.081 

25 53 42.027+1.9488 22 50 0.864+0.064 19 44 0.1340+0.013 

25 54 43.256(Intcrp.) 22 51 0.061+0.004 19 45 0.0196+0.002 

25 55 30.763+1.1455 22 52 0.005+0.001 18 35 0.0350+0.004 

25 56 0.280+0.0114 21 41 0.013+0.002 18 36 0.9479+0.107 

24 47 0.045+0.0041 21 42 0.581+0.061 18 37 3.8845+0.434 

24 48 0.792+0.0672 21 43 5.371+0.551 18 38 7.5135+0.827 

24 49 5.976+0.4774 21 44 13.897+1.375 18 39 4.7387+0.567 

24 50 23.198+1.7175 21 45 15.137+1.632 18 40 1.5857+0.188 

24 51 32.541+2.0109 21 46 8.244+0.918 18 41 0.3555+0.038 

24 52 39.839+3.6168 21 47 2.406+0.215 18 42 0.0639+0.006 

24 53 9.190+0.4258 21 48 0.397+0.034 18 43 0.0082+0.001 

24 54 2.733+0.1095 21 49 0.054+0.004 17 33 0.0154+0.002 

24 55 0.021+0.0010 21 50 0.002+0.001 17 34 0.4418+0.051 

23 45 0.032+0.0031 20 39 0.038+0.004 17 35 3.1153+0.357 

23 46 0.614+0.0568 20 40 1.148+0.123 17 36 4.0226+0.496 

23 47 5.827+0.5195 20 41 6.151+0.657 17 37 3.5159+0.430 



23 48 18.302+1.5607 20 42 12.461+1.296 17 38 0.7868+0.087 



500 A MeV 



Z A a(mb) 


Z 


A 


a(mb) 


17 39 188+0 0206 


13 30 0.083+0.010 


17 40 028±0 0031 


13 31 0.009+0.001 


1 fi 31 031 +0 0041 


12 


23 0.009+0.002 


16 32 822±0 0068 


12 


24 0.354+0.044 


16 33 2 686±0 3109 


12 25 0.748+0.093 


16 34 4 570±0 5764 


12 


26 


1.359+0.188 


16 35 1 862±0 2337 


12 


27 0.206+0.028 


16 36 499+0 0566 


12 


28 0.036+0.004 


16 37 070±0 0082 


11 


22 0.113+0.018 


16 38 012±0 0010 


11 


23 0.607+0.085 


15 29 015±0 0021 


11 


24 0.570+0.077 


15 30 378±0 0453 


11 25 0.199+0.025 


15 31 1 780±0 2119 


11 


26 0.024+0.003 


15 32 2.129±0.2726 


10 


22 


0.487+0.080 


15 33 1 081 +0 1 368 

1U iJ'J l.UOJ — LU . 1UUU 


10 


23 0.125+0.017 


15 34 219±0 0257 


10 24 0.014+0.002 


1 5 35 037+0 0041 








0.000+0.000 


1 5 36 003+0 001 








0.000+0.000 


14 27 026+0 0041 








0.000+0.000 


14 28 819+0 0977 








0.000+0.000 


14 29 1 500+0 1810 








0.000+0.000 


14 30 3 147+0 4103 








0.000+0.000 


14 31 0.614+0.0792 








0.000+0.000 


14 32 0.167+0.0195 








0.000+0.000 


14 33 0.013+0.0021 








0.000+0.000 


13 25 0.008+0.0010 








0.000+0.000 


13 26 0.254+0.0319 








0.000+0.000 


13 27 1.110+0.1459 








0.000+0.000 


13 28 1.311+0.1716 








0.000+0.000 


13 29 0.383+0.0504 








0.000+0.000 





300 A MeV 






a R (mb) = 701 ±56 




Z A 


rr ( Til hi 

Wilt V \J 1 


Z 


A 


rrim hi 


Z A rr(mb) 


zY 04 


n nnoj_n c\c\a~\ 

u.uyz±u.uu4i 


23 


50 


oo oin_i_o a 1 1 


zU 44 0./oO/±U.OoO 


Zl 00 


n m q_i_a HQ/ii 
U.yio±U.Uo41 


23 


51 


a nQK_i_n /ion 


on a c n oofinj-H nn/i 

zu 4o u.yoou±u.uy4 


Z7 00 
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23 


52 


i m a _i_n hkt 
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Of: k i 
ZD 01 


n noi _i_n nni n 


23 


53 


n i QQ_i_n (\(\r 
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zU 4/ U.UiyD±U.UUz 


Ofi c;o 
ZD Oz 


n c: k /i _i_n noin 
U.004±U.UolU 


23 54 


n nni _i_n nni 


i n Q7 n m i o_Ln nni 
iy of U.Ullo±U.UUl 


of? K O 

ZD Oo 


D.oU0±U.z40z 


22 


43 


n nQn_i_n nn/i 


i n qq n ot;nc_Ln non 

iy oo u.zDUo±u.uzy 


Of? K /I 

ZD 04 


oo.Uo4±l.ozyi 


22 


44 


n nioxn nno 


in on o ttoi _i_n ono 
iy oy z.MZlzhU.oUz 


26 55 
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22 


45 


^ a i o_i_n fin 
0.41z±U.0lz 


in /in k 1 Q/in i n £^o 
iy 4U 0.1o4y±U.00o 


25 49 


0.020±0.0021 
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zl.000±Z.UUU 


i n /i i /i ooooj_n nf^n 

iy 4i 4.oozo±u.yoy 


25 50 
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22 


47 
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iy 4z i.yyi/zhu.zoi 


25 51 
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zU.U ( o±z.Uz4 


1 n a o n cf;7n-Ln n^n 

iy 4o u.oo/u±u.uoy 


25 52 
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o. f OU±U.oUU 
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iy 44 U.Uo /DzhU.UUo 


25 53 
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22 


50 
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in a c nniO/i-i_n nni 
iy 40 U.Ulz4d_U.UUl 
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22 


51 
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lo oo U.UloOzhU.UUz 


25 55 


37.652±1.4019 


22 


52 
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U.UU4±U.UU1 
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lo oO U.40lO±U.Uoo 


25 56 
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21 


41 
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21 
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21 
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45 
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lo 41 U.zloo±U.Uzo 


24 51 


38.291±2.3654 


21 


46 


e.lof ±1. (54 


i o /io n not; i _i_n nn a 
lo 4z U.Uo0l±U.UU4 


24 52 


33.797±7.3473 


21 


47 
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24 53 


10.867±0.5043 


21 


48 


n Qc;o_i_n noi 
U.363±U.031 


1 t oo n nntio_Ln nni 
1 { oo U.UU0Z±U.UU1 


24 54 


2.847±0.1147 


21 49 


0.046±0.004 


17 34 0.1596±0.019 


24 55 


0.008±0.0010 


21 


50 


0.002±0.001 


17 35 1.4263±0.164 


23 45 


0.028±0.0031 


20 39 


0.020±0.002 


17 36 2.2440±0.255 


23 46 


0.640±0.0599 


20 40 


0.798±0.086 


17 37 1.6261±0.201 


23 47 


6.500±0.5795 


20 41 


4.648±0.497 


17 38 0.3388±0.038 


23 48 


21.757±1.8561 


20 42 


9.878±1.029 


17 39 0.0762±0.008 


23 49 


33.178±2.6514 


20 43 


9.838±1.836 


17 40 0.0134±0.001 



300 A MeV 

Z A a(mb) Z A a(mb) 
17 41 0.002±0.0010 14 31 0.386±0.051 
16 31 0.006±0.0010 14 32 0.047±0.006 
16 32 0.271±0.0319 14 33 0.005±0.001 
16 33 1.199±0.1390 13 26 0.059±0.007 
16 34 1.888±0.2213 13 27 0.438±0.059 
16 35 0.827±0.1050 13 28 0.491±0.067 
16 36 0.239±0.0268 13 29 0.334±0.044 
16 37 0.027±0.0031 13 30 0.098±0.013 
16 38 0.005±0.0010 13 31 0.004±0.001 
15 29 0.007±0.0010 12 24 0.023±0.004 
15 30 0.081±0.0103 12 25 0.288±0.037 
15 31 0.670±0.0802 12 26 0.399±0.056 
15 32 1.216±0.1574 12 27 0.120±0.016 
15 33 0.443±0.0576 12 28 0.004±0.001 
15 34 0.062±0.0072 11 23 0.280±0.040 
15 35 0.010±0.0010 11 24 0.319±0.041 
14 28 0.165±0.0195 11 25 0.078±0.010 
14 29 0.571±0.0689 10 22 0.316±0.053 



14 30 1.146±0.1522 10 23 0.041±0.006 



40 



VII. APPENDIX B: LONGITUDINAL VELOCITY OF RESIDUAL NUCLEI (MEAN AND R.M.S. 

VALUE) IN THE IRON BEAM SYSTEM AT REST. 







1500 A MeV 






1= Interpolated 


Z A Mean velocity(cm/: 


ns) R.M.S. (cm/ns) 


Z A Mean velocity(cm/: 


ns) R.M.S. (cm/ns) 


27 54 


-0.079±0.020 


0.127 


22 47 


-0.098(1) 


0.198 


27 55 


-0.058±0.015 


0.110 


22 48 


-0.089(1) 


0.176 


27 56 


-0.033±0.008 


0.104 


22 49 


-0.086±0.017 


0.172 


26 51 


-0.093±0.023 


0.141 


22 50 


-0.074±0.015 


0.161 


26 52 


-0.056±0.014 


0.132 


22 51 


-0.075±0.015 


0.165 


26 53 


-0.035±0.009 


0.114 


22 52 


-0.056±0.017 


0.138 


26 54 


-0.031±0.008 


0.084 


21 41 


-0.136±0.034 


0.275 


26 55 


-0.034±0.008 


0.072 


21 42 


-0.123±0.025 


0.273 


25 49 


-0.089±0.022 


0.177 


21 43 


-0.113±0.023 


0.251 


25 50 


-0.078±0.019 


0.168 


21 44 


-0.118±0.024 


0.237 


25 51 


-0.052±0.013 


0.150 


21 45 


-0.106(1) 


0.225 


25 52 


-0.053±0.013 


0.125 


21 46 


-0.104(1) 


0.208 


25 53 


-0.061±0.015 


0.105 


21 47 


-0.098±0.020 


0.199 


25 55 


-0.032±0.008 


0.074 


21 48 


-0.089±0.018 


0.192 


24 46 


-0.112±0.028 


0.201 


21 49 


-0.080±0.016 


0.172 


24 47 


-0.104±0.026 


0.191 


21 50 


-0.082±0.001 


0.157 


24 48 


-0.084±0.021 


0.189 


20 39 


-0.160±0.040 


0.293 


24 49 


-0.069±0.017 


0.181 


20 40 


-0.131±0.026 


0.288 


24 50 


-0.070±0.017 


0.167 


20 41 


-0.129±0.026 


0.274 


24 51 


-0.056±0.014 


0.160 


20 42 


-0.130±0.026 


0.260 


24 52 


-0.529(1) 


0.140 


20 43 


-0.132(1) 


0.259 


24 53 


-0.045±0.011 


0.117 


20 44 


-0.135(1) 


0.235 


24 54 


-0.023±0.006 


0.092 


20 45 


-0.116±0.023 


0.233 


24 55 


-0.040±0.022 


0.093 


20 46 


-0.105±0.021 


0.216 


23 45 


-0.118±0.030 


0.228 


20 47 


-0.102±0.020 


0.204 


23 46 


-0.095±0.019 


0.224 


20 48 


-0.097±0.019 


0.197 


23 47 


-0.084±0.017 


0.204 


19 37 


-0.167±0.033 


0.330 


23 48 


-0.088±0.018 


0.183 


19 38 


-0.145±0.022 


0.278 


23 49 


-0.081(1) 


0.168 


19 39 


-0.144±0.022 


0.297 


23 50 


-0.074(1) 


0.156 


19 40 


-0.145±0.022 


0.284 


23 51 


-0.064±0.013 


0.145 


19 41 


-0.137(1) 


0.284 


23 52 


-0.060±0.012 


0.121 


19 42 


-0.135(1) 


0.255 


23 53 


-0.061±0.015 


0.121 


19 43 


-0.121±0.018 


0.254 


22 43 


-0.123±0.016 


0.244 


19 44 


-0.119±0.018 


0.240 


22 44 


-0.107±0.018 


0.232 


19 45 


-0.112±0.022 


0.231 


22 45 


-0.097±0.019 


0.227 


19 46 


-0.122±0.025 


0.217 


22 46 


-0.103±0.021 


0.209 


18 35 


-0.157±0.024 


0.360 



1500 A MeV 1= Interpolated 



Z A Mean velocity (cm/: 


as) R.M.S. (cm/ns) 


Z 


A Mean velocity (cm/: 


ns) R.M.S. (cm/ns) 


18 36 


-0.164±0.025 


0.342 


13 


26 


-0.242±0.048 


0.510 


18 37 


-0.156±0.023 


0.322 


13 27 


-0.232±0.046 


0.467 


18 38 


-0.158±0.024 


0.291 


13 


28 


-0.221±0.044 


0.449 


18 39 


-0.148(1) 


0.285 


13 


29 


-0.220(1) 


0.448 


18 40 


-0.132(1) 


0.273 


13 30 


-0.206±0.041 


0.458 


18 41 


-0.136±0.020 


0.284 


13 


31 


-0.182±0.036 


0.445 


18 42 


-0.128±0.019 


0.261 


12 


24 


-0.225±0.045 


0.496 


18 43 


-0.128±0.019 


0.263 


12 


25 


-0.207±0.041 


0.487 


17 33 


-0.136±0.030 


0.420 


12 


26 


-0.20 (I) 


0.485 


17 34 


-0.177±0.026 


0.384 


12 


27 


-0.198±0.040 


0.465 


17 35 


-0.169±0.025 


0.345 


12 


28 


-0.218±0.054 


0.467 


17 36 


-0.158±0.024 


0.331 


11 


22 


-0.223±0.045 


0.561 


17 37 


-0.156(1) 


0.320 


11 


23 


-0.220±0.044 


0.550 


17 38 


-0.153±0.023 


0.308 


11 


24 


-0.218(1) 


0.533 


17 39 


-0.151±0.023 


0.294 


11 


25 


-0.214±0.043 


0.513 


17 40 


-0.146±0.022 


0.303 


11 


26 


-0.216±0.043 


0.504 


17 41 


-0.125±0.019 


0.294 


10 20 


-0.266±0.080 


0.630 


16 31 


-0.163±0.031 


0.412 


10 


21 


-0.260(1) 


0.599 


16 32 


-0.194±0.029 


0.397 


10 


22 


-0.250(1) 


0.583 


16 33 


-0.183±0.027 


0.375 


10 


23 


-0.237±0.071 


0.572 


16 34 


-0.186(1) 


0.355 


10 


24 


-0.153±0.095 


0.569 


16 35 


-0.180(1) 


0.351 


9 


18 


-0.253±0.076 


0.680 


16 36 


-0.183±0.027 


0.339 


9 


19 


-0.260±0.078 


0.650 


16 37 


-0.176±0.026 


0.345 


9 


20 


-0.260±0.078 


0.656 


16 38 


-0.165±0.025 


0.322 


9 


21 


-0.225±0.067 


0.641 


15 30 


-0.213±0.032 


0.439 


8 


16 


-0.240±0.072 


0.743 


15 31 


-0.201±0.030 


0.401 


8 


17 


-0.230±0.069 


0.715 


15 32 


-0.198(1) 


0.391 


8 


18 


-0.250±0.075 


0.703 


15 33 


-0.197(1) 


0.379 












n i QR+n D9Q 

u.iyuxu.uiy 


0.379 










15 35 


-0.187±0.028 


0.364 










15 36 


-0.162±0.024 


0.348 










14 28 


-0.221±0.033 


0.447 










14 29 


-0.218±0.033 


0.436 










14 30 


-0.219(1) 


0.412 










14 31 


-0.210(1) 


0.412 










14 32 


-0.202±0.030 


0.405 










14 33 


-0.202±0.030 


0.402 















1000 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


nsl RMS frm/ns'l 


26 51 


-0.086±0.022 


0.144 


22 48 


-0.095(1) 


0.176 


26 52 


-0.075±0.019 


a i i t 
0.11 1 


22 49 


-0.083±0.017 


0.162 


26 53 


-0.056±0.014 


0.099 


22 50 


-0.073±0.015 


0.146 


26 54 


-0.041±0.010 


O.Ooo 


22 51 


-0.074±0.015 


0.140 


26 55 


-0.031±0.008 


a acc 
0.05b 


22 52 


-0.064±0.013 


0.126 


25 48 


-0.097±0.024 


a i oo 
0.182 


21 40 


-0.166±0.033 


0.2(4 


25 49 


-0.097±0.024 


a i to 
0.1 (8 


21 41 


-0.157±0.031 


0.260 


25 50 


-0.099±0.025 


0.158 


21 42 


-0.147±0.029 


0.2o» 


25 51 


-0.078±0.019 


a i oo 
0.138 


21 43 


-0.138±0.028 


A OOA 


25 52 


-0.068±0.017 


a i i o 
0.118 


21 44 


-0.127±0.025 


A 001 


25 53 


-0.053±0.013 


a aaq 
0.098 


21 45 


-0.121(1) 


A OO 1 


25 55 


-0.049±0.015 


a nco 
0.008 


21 46 


-0.110(1) 


A OAO 

O.zUo 


24 46 


-0.114±0.028 


a oak 
0.2U5 


21 47 


-0.099±0.020 


A 1 A1 

o.iyi 


24 47 


-0.106±0.027 


0.199 


21 48 


-0.094±0.019 


A 1 O 1 

0.181 


24 48 


-0.109±0.027 


0.1 to 


21 49 


-0.085±0.017 


A 1 dQ 

0.168 


24 49 


-0.092±0.023 


n i tip. 
O.lDD 


21 50 


-0.089±0.022 


A 1 C O 

0.152 


24 50 


-0.085±0.021 


0.143 


20 38 


-0.134±0.034 


0.148 


24 51 


-0.074±0.019 


0.130 


20 39 


-0.167±0.033 


0.293 


24 52 


-0.062(1) 


0.112 


20 40 


-0.159±0.032 


0.272 


24 53 


-0.047±0.012 


0.104 


20 41 


-0.151±0.030 


0.262 


24 54 


-0.036±0.009 


0.085 


20 42 


-0.144±0.029 


0.000 


23 44 


-0.124±0.031 


0.000 


20 43 


-0.135(1) 


0.245 


23 45 


-0.127±0.025 


0.235 


20 44 


-0.124(1) 


0.229 


23 46 


-0.124±0.025 


0.208 


20 45 


-0.119±0.024 


0.218 


23 47 


-0.108±0.022 


0.193 


20 46 


-0.107±0.021 


0.204 


23 48 


-0.097±0.019 


0.174 


20 47 


-0.107±0.021 


0.193 


23 49 


-0.090±0.018 


0.158 


20 48 


-0.095±0.019 


0.168 


23 50 


-0.079(1) 


0.150 


19 36 


-0.149±0.030 


0.000 


23 51 


-0.067±0.013 


0.133 


19 37 


-0.158±0.032 


A O 1 1 

0.311 


23 52 


-0.056±0.011 


0.122 


19 38 


-0.167±0.033 


0.305 


23 53 


-0.049±0.010 


0.110 


19 39 


-0.172±0.034 


0.282 


23 54 


-0.031±0.009 


0.102 


19 40 


-0.152±0.030 


0.275 


22 43 


-0.139±0.028 


0.240 


19 41 


-0.151(1) 


0.268 


22 44 


-0.136±0.027 


0.227 


19 42 


-0.143(1) 


0.256 


22 45 


-0.123±0.025 


0.216 


19 43 


-0.132±0.020 


0.243 


22 46 


-0.112±0.022 


0.195 


19 44 


-0.120±0.018 


0.232 



22 47 



-0.104(1) 



0.195 19 45 -0.114±0.017 



0.213 







1000 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


ns*) RMS frm/ns'l 


18 35 


-0.181±0.036 


O.oo9zo 


14 30 


-0.207±0.031 


0.405 


18 36 


-0.184±0.028 


U.ol / ZD 


14 31 


-0.203±0.030 


0.390 


18 37 


-0.179±0.027 


r» oi 1 OP 

U.ollzo 


14 32 


-0.192±0.029 


0.382 


18 38 


-0.164±0.025 


U.ZyUzO 


14 33 


-0.191±0.029 


0.361 


18 39 


-0.158(1) 


U.zozzo 


13 24 


-0.169±0.042 


n £;o a 
0.624 


18 40 


-0.142±0.021 


U.ZDDZO 


13 25 


-0.188±0.038 


0.01 l 


18 41 


-0.143±0.022 


U.ZOZZO 


13 26 


-0.210±0.042 


n A O A 

0.494 


18 42 


-0.141±0.021 


U.ZOOZO 


13 27 


-0.207±0.041 


0.456 


18 43 


-0.130±0.020 


U.zoUzo 


13 28 


-0.212(1) 


n A A i 
0.441 


17 32 


0.000±0.000 


n o£! cok 
U.OOOZO 


13 29 


-0.208(1) 


n /loo 
0.42.5 


17 33 


-0.186±0.037 


U.OOZZO 


13 30 


-0.206±0.041 


n A T7 

0.42 I 


17 34 


-0.189±0.028 


U.OOZZO 


13 31 


-0.203(1) 


0.393 


17 35 


-0.181±0.027 


U.oo4z4 


12 22 


-0.171±0.051 


0.069 


17 36 


-0.178±0.027 


n o i i o a 
U.ollz4 


12 23 


-0.208±0.052 


n nnn 
0.000 


17 37 


-0.173(1) 


n onoo/i 
U.oUoz4 


12 24 


-0.219±0.044 


0.021 


17 38 


-0.163±0.024 


n on aha 
U.zy4z4 


12 25 


-0.211±0.042 


n a da 
0.496 


17 39 


-0.158±0.024 


n ono o /i 

U.z9oz4 


12 26 


-0.210(1) 


0.466 


17 40 


-0.152±0.023 


C\ 070 O /I 

(J.z7oz4 


12 27 


-0.207(1) 


0.451 


16 31 


-0.193±0.039 


0.39724 


12 28 


-0.208±0.042 


0.455 


16 32 


-0.201±0.030 


C\ O TOO /I 

U.o7zz4 


11 20 


0.000±0.000 


0.000 


16 33 


-0.200±0.030 


0.36024 


11 21 


-0.184±0.055 


0.611 


16 34 


-0.189(1) 


0.33523 


11 22 


-0.209±0.042 


0.598 


16 35 


-0.183(1) 


0.33423 


11 23 


-0.208±0.042 


0.547 


16 36 


-0.172±0.026 


n oi coo 

0.31623 


11 24 


-0.203±0.041 


0.511 


16 37 


-0.167±0.025 


U.ol4zo 


11 25 


-0.197±0.039 


0.503 


16 38 


-0.157±0.024 


C\ O 1 1 o o 

U.ollzo 


11 26 


-0.191±0.038 


0.493 


15 28 


-0.175±0.044 


0.00023 


10 19 


-0.185±0.055 


0.652 


15 29 


-0.210±0.031 


r» di coo 

U.41ozo 


10 20 


-0.223±0.067 


0.627 


15 30 


-0.207±0.031 


n /i i noo 
U.41Uzo 


10 21 


-0.210(1) 


0.00( 


15 31 


-0.203±0.030 


n ooooo 

u.oyzzo 


10 22 


-0.209(1) 


0.000 


15 32 


-0.197(1) 


u.o / izo 


10 23 


-0.257±0.077 


U.040 


15 33 


-0.196(1) 


0.35923 


9 18 


-0.217±0.065 


0.687 


15 34 


-0.190±0.028 


0.35222 


9 19 


-0.230(1) 


0.639 


15 35 


-0.184±0.028 


0.34222 


9 20 


-0.220(1) 


0.624 


14 27 


-0.194±0.049 


0.47522 


9 21 


-0.191±0.057 


0.647 


14 28 


-0.215±0.032 


0.43222 


8 17 


-0.250(1) 


0.692 


14 29 


-0.208±0.031 


0.41822 


8 18 


-0.240(1) 


0.680 







750 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


nsl RMS fcm/nsl 


27 54 


-0.060±0.015 


0.096 


22 47 


-0.122(1) 


0.199 


27 55 


-0.058±0.015 


n hto 
0.0/3 


22 48 


-0.106(1) 


n i *7 A 
0.1 (4 


27 56 


-0.049±0.012 


0.058 


22 49 


-0.101±0.020 


0.169 


26 51 


-0.083±0.021 


0.130 


22 50 


-0.084±0.017 


n i Kn 
O.loO 


26 52 


-0.053±0.013 


0.118 


22 51 


-0.089±0.018 


n i A A 
U.144 


26 53 


-0.050±0.013 


n i no 
0.102 


22 52 


-0.083±0.025 


n 1 on 
0.130 


26 54 


-0.043±0.011 


O.Oos 


21 41 


-0.111±0.033 


n o a c 
0.240 


26 55 


-0.036±0.009 


0.003 


21 42 


-0.136±0.027 


n no c 
0.2JO 


25 49 


-0.093±0.023 


0.1 /9 


21 43 


-0.143±0.029 


n ion 
U.2o2 


25 50 


-0.087±0.022 


f\ 1 C A 

0.154 


21 44 


-0.153±0.031 


n oon 
0.2o0 


25 51 


-0.075±0.019 


n i q A 
0.134 


21 45 


-0.145(1) 


n ooc; 
U.ZZO 


25 52 


-0.076±0.019 


0.11 1 


21 46 


-0.131(1) 


O.zUo 


25 53 


-0.063±0.016 


n i n a 
0.104 


21 47 


-0.123±0.025 


n i no 
0.193 


25 55 


-0.022±0.005 


0.00/ 


21 48 


-0.109±0.022 


n i oo 
0.183 


24 47 


-0.107±0.027 


n i nn 
0.190 


21 49 


-0.105±0.021 


n i to 
0.1 /3 


24 48 


-0.093±0.023 


0.181 


21 50 


-0.106±0.027 


n i c t 
0.10 / 


24 49 


-0.094±0.023 


0.163 


20 40 


-0.162±0.049 


0.266 


24 50 


-0.094±0.024 


0.145 


20 41 


-0.159±0.032 


0.250 


24 51 


-0.089±0.022 


0.133 


20 42 


-0.167±0.033 


0.249 


24 52 


-0.072(1) 


0.115 


20 43 


-0.158(1) 


0.250 


24 53 


-0.053±0.013 


0.107 


20 44 


-0.150(1) 


0.231 


24 54 


-0.041±0.010 


0.091 


20 45 


-0.145±0.029 


0.223 


24 55 


-0.030±0.015 


0.092 


20 46 


-0.128±0.026 


0.209 


23 45 


-0.136±0.041 


0.211 


20 47 


-0.121±0.024 


0.194 


23 46 


-0.117±0.023 


0.210 


20 48 


-0.122±0.037 


0.186 


23 47 


-0.111±0.022 


0.188 


19 38 


-0.199±0.040 


0.280 


23 48 


-0.117±0.023 


0.178 


19 39 


-0.189±0.028 


0.277 


23 49 


-0.109±0.022 


0.169 


19 40 


-0.184±0.028 


0.274 


ZO JU 


n nsfi/n 

u-UOU^l 1 


n iro 

0.153 


19 41 


n i 7Qfn 


U.ZOZ 


23 51 


-0.073±0.015 


0.13/ 


19 42 


-0.164(1) 


0.ZO9 


23 52 


-0.072±0.014 


0.126 


19 43 


-0.158±0.024 


0.250 


23 53 


-0.060±0.012 


0.111 


19 44 


-0.148±0.022 


0.236 


23 54 


-0.119±0.036 


0.102 


19 45 


-0.136±0.020 


0.218 


22 43 


-0.130±0.039 


0.231 


19 46 


-0.142±0.028 


0.212 


22 44 


-0.123±0.025 


0.233 


18 36 


-0.156±0.035 


0.306 


22 45 


-0.127±0.025 


0.209 


18 37 


-0.178±0.036 


0.292 


22 46 


-0.133±0.027 


0.202 


18 38 


-0.184±0.028 


0.293 







750 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


nsl RMS frm/ns'l 


18 39 


-0.169(1) 


0.285 


13 31 


-0.228±0.057 


0.399 


18 40 


-0.181±0.027 


n 0£?n 
0.2D0 


12 24 


-0.140±0.072 


n c 1 1 
O.oll 


18 41 


-0.179±0.027 


0.250 


12 25 


-0.226±0.056 


0.470 


18 42 


-0.161±0.024 


n o a n 
0.249 


12 26 


-0.230(1) 


0.4D8 


18 43 


-0.162±0.032 


n o a o 
0.242 


12 27 


-0.259±0.065 


0.458 


17 34 


-0.152±0.045 


n o k o 
0.358 


12 28 


-0.214±0.064 


n a An 
0.449 


17 35 


-0.215±0.043 


nolo 
0.312 


11 22 


-0.123±0.071 


n coi 
0.031 


17 36 


-0.201(1) 


0.316 


11 23 


-0.238±0.072 


n coi 
0.021 


17 37 


-0.190(1) 


n one 
0.305 


11 24 


-0.220(1) 


n k i n 
0.010 


17 38 


-0.183±0.027 


n onn 
0.290 


11 25 


-0.253±0.076 


O.olb 


17 39 


-0.172±0.026 


n oon 
0.289 


11 26 


-0.209±0.063 


n /i o /i 
0.484 


17 40 


-0.167±0.025 


n otk 
0.2/0 


10 21 


-0.246(1) 


n /i 7n 
0.4(0 


16 32 


-0.188±0.042 


n 0£? A 

0.3b4 


10 22 


-0.240(1) 


n Knn 
O.oOO 


16 33 


-0.202±0.040 


n o a o 
0.348 


10 23 


-0.299±0.090 


n ko 1 
0.031 


16 34 


-0.202(1) 


n oon 
0.339 


10 24 


-0.193±0.058 


n c 1 o 
0.012 


16 35 


-0.195(1) 


n ooc; 
0.32b 


9 19 


-0.293(1) 


n k o k 
0.080 


16 36 


-0.202±0.030 


0.316 


9 20 


-0.280(1) 


0.569 


16 37 


-0.190±0.028 


0.314 


8 17 


-0.300(1) 


0.629 


16 38 


-0.173±0.035 


0.285 


8 18 


-0.310(1) 


0.606 


15 30 


-0.153±0.051 


0.402 








15 31 


-0.211±0.042 


0.377 








15 32 


-0.210(1) 


0.360 








15 33 


-0.202(1) 


0.356 








15 34 


-0.227±0.045 


0.358 








15 35 


-0.187±0.037 


0.330 








15 36 


-0.201±0.040 


0.311 








14 29 


-0.212±0.053 


0.411 








14 30 


-0.223±0.045 


0.403 








14 31 


n OOA 1 Ci C\ A A 

-0.220±0.044 


n ooo 
0.382 








14 32 


-0.219±0.044 


n ooo 
0.382 








14 33 


-0.195±0.039 


fl OKQ 








14 34 


-0.198±0.049 


0.347 








13 26 


-0.240±0.067 


0.431 








13 27 


-0.223±0.056 


0.437 








13 28 


-0.228(1) 


0.437 








13 29 


-0.224(1) 


0.396 








13 30 


-0.246±0.062 


0.407 













500 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


ns*) RMS frm/ns'l 


27 53 


-0.097±0.013 


0.106 


22 47 


-0.156(1) 


0.193 


27 54 


-0.088±0.012 


n nnn 
U.090 


22 48 


-0.149(1) 


n i T/i 
0.1 (4 


27 55 


-0.073±0.012 


0.067 


22 49 


-0.106±0.021 


0.169 


27 56 


-0.074±0.015 


U.058 


22 50 


-0.092±0.018 


n 1 ko 
0.153 


26 51 


-0.076±0.019 


U.14D 


22 51 


-0.090±0.018 


n 1 (In 

U.149 


26 52 


-0.066±0.016 


U.121 


22 52 


-0.081±0.016 


0.131 


26 53 


-0.049±0.012 


0.100 


21 41 


-0.195±0.039 


n one 
0.305 


26 54 


-0.034±0.008 


0.0/0 


21 42 


-0.188±0.038 


0.2o 1 


26 55 


-0.027±0.007 


O.Ooe 


21 43 


-0.176±0.035 


n o A i 
U.241 


25 49 


-0.122±0.030 


0.191 


21 44 


-0.162±0.032 


n oi7 
0.22 1 


25 50 


-0.103±0.026 


n i £?o 


21 45 


-0.179(1) 


n oo i 
U.zzl 


25 51 


-0.081±0.020 


n i A i 
(J. 141 


21 46 


-0.166(1) 


n ono 
O.zllo 


25 52 


-0.066±0.016 


U.121 


21 47 


-0.1880±0.026 


n i nT 
0.19 1 


25 53 


-0.057±0.014 


n i no 
U.108 


21 48 


-0.120±0.024 


n i 

0.186 


25 55 


-0.024±0.006 


0.0/4 


21 49 


-0.088±0.018 


n i T/i 
0.1 (4 


25 56 


-0.160±0.080 


U.ObO 


21 50 


-0.095±0.024 


rt i An 
U.142 


24 47 


-0.137±0.034 


0.207 


20 39 


-0.220±0.055 


0.300 


24 48 


-0.124±0.031 


0.183 


20 40 


-0.192±0.038 


0.273 


24 49 


-0.109±0.027 


0.170 


20 41 


-0.186±0.037 


0.266 


24 50 


-0.090±0.023 


0.148 


20 42 


-0.182±0.036 


0.243 


24 51 


-0.092±0.023 


0.139 


20 43 


-0.175(1) 


0.236 


24 52 


-0.068(1) 


0.118 


20 44 


-0.160(1) 


0.229 


24 53 


-0.049±0.012 


0.109 


20 45 


-0.149±0.030 


0.223 


24 54 


-0.041±0.010 


0.095 


20 46 


-0.134±0.027 


0.208 


23 45 


-0.171±0.034 


0.233 


20 47 


-0.104±0.026 


0.196 


23 46 


-0.150±0.030 


0.210 


19 37 


-0.225±0.045 


0.322 


23 47 


-0.131±0.026 


0.196 


19 38 


-0.216±0.032 


0.301 


23 48 


-0.119±0.024 


0.178 


19 39 


-0.199±0.030 


0.286 


23 49 


-0.117±0.023 


0.109 


19 40 


a i no i n nin 

-U.193±0.029 


n o^c 
0.2ob 


23 50 


-0.094(1) 


n i co 
U.152 


19 41 


-0.183(1) 


n oco 
O.ZOO 


23 51 


-0.077±0.015 


n i /in 


19 42 


-0.175(1) 


U.zOO 


23 52 


-0.071±0.014 


0.129 


19 43 


-0.165±0.025 


0.248 


23 53 


-0.064±0.013 


0.118 


19 44 


-0.156±0.023 


0.238 


22 43 


-0.185±0.037 


0.258 


19 45 


-0.120±0.024 


0.207 


22 44 


-0.164±0.033 


0.230 


19 46 


O.OOOiO.OOO 


0.000 


22 45 


-0.156±0.031 


0.218 


18 35 


-0.243±0.049 


0.342 


22 46 


-0.141±0.028 


0.201 


18 36 


-0.221±0.033 


0.320 







500 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


nsl RMS fcm/nsl 


18 37 


-0.208±0.031 


0.313 


13 26 


-0.258±0.052 


0.522 


18 38 


-0.205±0.031 


n onK 
0.295 


13 27 


-0.234±0.047 


n a p.q 
0.456 


18 39 


-0.200(1) 


0.283 


13 28 


-0.230(1) 


0.421 


18 40 


-0.196(1) 


n oc i 
0.251 


13 29 


-0.228(1) 


n A n*7 
0.40/ 


18 41 


-0.190±0.028 


0.2/4 


13 30 


-0.223±0.045 


n a 7*0 
0.4/2 


18 42 


-0.151±0.023 


0.253 


12 24 


-0.287±0.057 


n c cn 
0.55U 


18 43 


-0.136±0.027 


A Ol O 

0.218 


12 25 


-0.268±0.054 


0.4/0 


17 33 


-0.236±0.047 


n oo i 
0.381 


12 26 


-0.260(1) 


n /i o k 
0.420 


17 34 


-0.225±0.034 


0.3b2 


12 27 


-0.256(1) 


n /i /i t 
0.44/ 


17 35 


-0.216±0.032 


n oo a 
0.334 


12 28 


-0.207±0.041 


n /loo 
0.482 


17 36 


-0.215(1) 


A OOO 

0.616 


11 22 


-0.297±0.074 


n cno 
0.093 


17 37 


-0.213(1) 


n one 


11 23 


-0.278±0.080 


n c ^n 
0.0 /O 


17 38 


-0.212±0.032 


n 007 
0.28/ 


11 24 


-0.265(1) 


n kt/o 
0.0 / 3 


17 39 


-0.205±0.031 


n ono 
0.302 


10 22 


-0.280(1) 


n koo 
0.028 


17 40 


-0.127±0.038 


n oon 
0.280 


10 23 


-0.241±0.100 


n cn7 
0.59 / 


16 31 


-0.286±0.086 


n oon 
0.389 


10 24 


-0.290±0.102 


n c 7"7 
0.0 / / 


16 32 


-0.240±0.036 


0.381 








16 33 


-0.236±0.034 


0.358 








16 34 


-0.231(1) 


0.331 








16 35 


-0.228(1) 


0.318 








16 36 


-0.223±0.033 


0.326 








16 37 


-0.209±0.031 


0.349 








16 38 


-0.128±0.038 


0.275 








15 29 


-0.257±0.051 


0.386 








15 30 


-0.234±0.035 


0.432 








15 31 


-0.232±0.035 


0.377 








15 32 


-0.224(1) 


0.385 








15 33 


-0.220(1) 


0.349 










n 91 R+n pur 


n ocn 
0.309 








15 35 


-0.262±0.058 


n OTn 
0.3/9 








14 27 


-0.303±0.061 


0.429 








14 28 


-0.266±0.040 


0.456 








14 29 


-0.251±0.038 


0.407 








14 30 


-0.246(1) 


0.405 








14 31 


-0.238(1) 


0.372 








14 32 


-0.225±0.034 


0.408 








14 33 


-0.213±0.043 


0.373 













300 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


nsl RMS fcm/nsl 


27 54 


-0.069±0.015 


0.084 


23 54 


-0.210(1) 


0.250 


27 55 


-0.060±0.012 


n nc c 
0.065 


22 43 


-0.213±0.043 


n oi n 
0.21U 


27 56 


-0.051±0.011 


0.058 


22 44 


-0.206±0.041 


0.224 


26 51 


-0.110±0.023 


n i of? 
0.136 


22 45 


-0.186±0.037 


n onn 
0.2U9 


26 52 


-0.082±0.021 


n i i c 
0.115 


22 46 


-0.178±0.036 


nine 


26 53 


-0.071±0.018 


n i nc 
0.106 


22 47 


-0.135(1) 


n i n /i 
U.194 


26 54 


-0.053±0.013 


n nn a 
0.094 


22 48 


-0.149(1) 


n i to 
0.1 11 


26 55 


-0.030±0.007 


n nTO 
0.0/3 


22 49 


-0.114±0.035 


n i a a 
0.164 


25 49 


-0.153±0.038 


n i 

0.162 


22 50 


-0.119±0.030 


n i c *7 
O.lo / 


25 50 


-0.128±0.032 


n i Kn 
0.150 


22 51 


-0.128±0.030 


n i A 'y 
0.14/ 


25 51 


-0.105±0.026 


n i o c 
0.135 


21 42 


-0.231±0.046 


n o ^ o 
0.248 


25 52 


-0.093±0.023 


n i on 
0.130 


21 43 


-0.204±0.041 


n 007 
0.2J / 


25 53 


-0.075±0.019 


n i n*7 
0.10/ 


21 44 


-0.196±0.039 


n ooc 


25 55 


-0.024±0.006 


n nTO 
0.0/2 


21 45 


-0.188(1) 


noil 


24 47 


-0.177±0.044 


n i oo 
0.182 


21 46 


-0.177(1) 


n om 
0.2U1 


24 48 


-0.155±0.039 


n i 

0.1 / / 


21 47 


-0.159±0.032 


n i 

0.186 


24 49 


-0.137±0.034 


0.162 


21 48 


-0.149±0.030 


0.186 


24 50 


-0.118±0.029 


0.154 


21 49 


-0.149±0.034 


0.172 


24 51 


-0.096±0.024 


0.139 


20 40 


-0.243±0.049 


0.264 


24 52 


-0.090(1) 


0.131 


20 41 


-0.224±0.045 


0.257 


24 53 


-0.076±0.019 


0.103 


20 42 


-0.206±0.041 


0.247 


24 54 


-0.063±0.016 


0.101 


20 43 


-0.207(1) 


0.238 


23 46 


-0.190±0.038 


0.201 


20 44 


-0.196(1) 


0.221 


23 47 


-0.162±0.032 


0.188 


20 45 


-0.189±0.038 


0.215 


23 48 


-0.148±0.030 


0.180 


20 46 


-0.164±0.033 


0.202 


23 49 


-0.125±0.025 


0.170 


20 47 


-0.156±0.031 


0.195 


23 50 


-0.120(1) 


0.155 


19 38 


-0.250±0.050 


0.285 


23 51 


-0.115±0.023 


0.138 


19 39 


-0.236±0.035 


0.282 


23 52 


-0.103±0.021 


0.135 


19 40 


-0.229±0.034 


0.269 



23 53 -0.096±0.019 



0.124 19 41 



-0.224(1) 



0.261 







300 A MeV 






1= Interpolated 


Z A Mean velocity (cm/: 


nsl RMS (cm /ns~l 


Z A Mean velocity (cm/: 


ns*) RMS frm/ns'l 


19 42 


-0.212(1) 


0.243 


15 35 


-0.262±0.039 


0.321 


19 43 


-0.210±0.031 


n o a o 
U.242 


14 28 


-0.308±0.062 


n A i o 
0.413 


19 44 


-0.184±0.037 


0.217 


14 29 


-0.305±0.046 


0.389 


19 45 


-0.209±0.042 


n om 
U.201 


14 30 


-0.300(1) 


n 070 
0.3(8 


18 36 


-0.261±0.039 


n oi A 
U.314 


14 31 


-0.295(1) 


nolo 
0.319 


18 37 


-0.255±0.038 


n qht 
U.29/ 


14 32 


-0.287±0.064 


n ono 
0.392 


18 38 


-0.200±0.042 


n oto 


13 26 


-0.339±0.102 


n A A T 

0.44/ 


18 39 


-0.236(1) 


n o^n 
0.2/9 


13 27 


-0.326(1) 


n /i nn 
0.400 


18 40 


-0.215±0.032 


n Of? a 
U.ZD4 


13 28 


-0.310(1) 


n oof; 
0.33b 


18 41 


-0.212±0.032 


ri 0*70 


13 29 


-0.308±0.082 


n of? c 
0.3OO 


18 42 


-0.212±0.042 


n ooo 
U.238 


12 25 


-0.325±0.098 


n oo c 
0.350 


17 34 


-0.272±0.054 


n oot 


12 26 


-0.318(1) 


n on /i 
0.394 


17 35 


-0.252±0.038 


n oo k 
U.335 


12 27 


-0.310±0.098 


n one 
0.390 


17 36 


-0.230(1) 


n ono 
0.292 


12 28 


-0.309±0.093 


n /in'? 
0.4U ( 


17 37 


-0.255(1) 


n one 
0.295 


11 23 


-0.326±0.098 


n oo /i 
0.384 


17 38 


-0.244±0.037 


n ooo 


11 24 


-0.310(1) 


n oc; e 
0.300 


17 39 


-0.188±0.043 


0.274 








17 40 


-0.243±0.036 


0.265 








16 32 


-0.298±0.060 


0.367 








16 33 


-0.266±0.040 


0.344 








16 34 


-0.253(1) 


0.333 








1U oo 




0.322 








16 36 


-0.280±0.065 


0.320 








16 37 


-0.255±0.038 


0.306 








16 38 


-0.223±0.045 


0.288 








15 30 


-0.322±0.064 


0.398 








15 31 


-0.295±0.044 


0.375 








15 32 


-0.290(1) 


0.370 








15 33 


-0.280(1) 


0.348 








15 34 


-0.272±0.041 


0.344 
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